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SUMMARY 


The  regular  (RR)  to  Mach  reflection  (MR)  transition  boundary  in 
pseudostationary  flow  in  air  was  investigated,  with  emphasis  on  the  influence 
of  the  shock-induced  boundary  layer.  Nearly  100  experiments  were  conducted  in 
the  region  of  the  RR-*-*41R  transition  line  over  a  range  of  incident-shock-wave 
Mach  number  1.1  <  M  <  6.5  by  conducting  experiments  in  the  UTIAS  10  x  18  cm 
Hypervelocity  Shock  Tube.  The  wedge  angles  used  were  42  ,  45  ,  47  ,  and  48  . 
Initial  pressures  were  kept  as  low  as  possible  to  maximize  viscous  effects  and 
ranged  from  0.4  kPa  to  100  kPa.  A  comparison  was  made  between  the 
boundary-layer  thickness,  and  the  deviation  of  the  transition  boundary  from 
inviscid  theory.  This  difference  is  often  called  the  'von  Neumann  paradox'. 

It  was  found  that  the  'von  Neumann  paradox'  was  due  to  viscous  effects. 
The  deviation  of  the  RR-*^MR  transition  line  from  the  'detachment'  criterion 
boundary  was  found  to  increase  with  a  drop  in  initial  pressure,  in  a  manner 
consistent  with  boundary-layer  theory.  The  effect  of  the  endwall  boundary 
layer  on  the  RR-**-MR  transition  line  was  more  pronounced  at  low  Mach  number 
(Ms  <  2),  and  a  model  is  proposed  to  explain  this  behavior. 

Initial  pressure  was  also  found  to  influence  the  height  of  the  Mach  stem 
in  MR.  Lower  initial  pressure  (and  hence  greater  viscous  effects)  reduced  the 
height  of  the  Mach  stem,  which  was  found  to  be  smaller  than  the  inviscid 
prediction  in  all  MR  experiments.  An  explanation  for  the  reduction  in  Mach 
stem  height  is  suggested,  but  the  cause  was  not  thoroughly  investigated. 
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Chapter  1 
INTRODUCTION 

When  energy  is  instantaneously  released  in  a  material,  an  explosion 
occurs  and  a  shock  wave  is  generated.  The  shock  wave,  which  causes  a  sudden 
jump  in  pressure,  temperature,  and  entropy,  travels  through  the  material  as  a 
result  of  molecular  collisions.  It  loses  strength  by  heating  the  material  and 
raising  its  pressure  as  it  races  outward  from  the  center  of  the  explosion,  and 
ultimately  weakens  to  a  sound  wave. 

If  a  planar  shock  wave  intersects  a  wedge  surface,  four  different 
reflection  patterns  can  occur,  depending  on  the  strength  of  the  shock  and  the 
wedge  angle.  The  patterns  are  regular  reflection  (RR),  single-Mach  reflection 
(SMR),  complex-Mach  reflection  (CMR),  and  double-Mach  reflection  (DMR)  (see 
Figs.  1  to  4,  from  Ref.  30).  The  latter  three  as  a  group  are  called  Mach 
reflection  (MR).  The  four  types  have  flow  fields  with  different  pressure, 
density  and  temperature  distributions  associated  with  them.  The  shock  wave 
from  an  above-ground  explosion  in  air  would  experience  several  of  the 
reflection  patterns  as  it  intersected  the  ground  on  its  travel  outward  [1,28]. 

There  are  many  examples  of  both  man-made  and  naturally  occurring 
shock-wave  phenomena  [2],  In  several  of  these,  the  process  of 
oblique-shock-wave  reflection  is  of  great  concern.  For  example,  accurate 
prediction  of  the  type  of  reflection  and  associated  pressures  is  necessary  for 
the  design  of  shock-resistant  structures.  Applications  range  from  chemical 
plants,  where  the  potential  for  major  explosions  always  exists,  to  mobile 
missile  launchers  which  must  be  able  to  withstand  a  nuclear  blast  [3], 


The  study  of  oblique-shock-wave  reflections  began  as  far  back  as  1878 
when  Ernst  Mach  [4]  did  research  on  shock  waves  and  observed  RR  and  MR.  In 
spite  of  this,  it  was  not  until  the  early  1940's  that  more  thorough  analytical 
and  experimental  work  was  done.  At  this  time  an  analytical  criterion  for  the 
transition  from  RR  to  MR  in  nonstationary  flow  (usually  referred  to  as  the 
'detachment'  criterion)  was  suggested  by  von  Neumann  [5,6]. 

Under  the  direction  of  Bleakney  [7]  a  major  study  of  nonstationary 
oblique-shock-wave  reflection  was  carried  out  at  Princeton  University.  Smith 
[8]  did  an  extensive  series  of  experiments  and  discovered  SMR  and  CMR.  He 
also  found  that  RR  persisted  beyond  the  boundary  predicted  by  the  'detachment' 
criterion.  This  persistence  is  often  referred  to  as  the  'von  Neumann 
paradox'.  Further  experimentation  by  White  [9]  confirmed  this  persistence. 
White  [9]  also  suggested  a  criterion  for  the  transition  from  SMR  to  CMR,  and 
discovered  DMR,  thus  completing  the  identification  of  the  four  major  types  of 
oblique-shock-wave  reflection.  The  most  current  analytical  transition 
boundaries  for  inviscid  frozen  air  are  shown  in  Fig.  5. 

Much  work  has  been  done  since  then  in  defining  the  regions  of  RR,  SMR, 
CMR,  and  DMR.  A  study  by  Kawamura  and  Saito  [10]  again  confirmed  the 
persistence  of  RR  into  the  predicted  MR  region  for  nonstationary  flows. 
Gvozdeva  et  al  [11]  and  Henderson  and  Lozzi  [12]  suggested  criteria  for  the 
transition  from  SMR  to  CMR,  and  CMR  to  DMR  but  they  were  not  substantiated. 
In  the  same  study,  Henderson  and  Lozzi  [12]  also  offered  the  'mechanical 
equilibrium'  criterion  for  transition  from  RR  to  MR  in  nonstationary  flows. 
Bazhenova  et  al  [13]  attempted  to  clarify  experimentally  the  regions  of  RR, 
SMR,  CMR,  and  DMR,  but  like  other  studies  their  range  of  Mach  numbers  and 
wedge  angles  was  limited.  Law  and  Glass  [14]  did  experiments  over  a  wider 


range,  and  later  Ben-Dor  and  Glass  [15]  established  the  generally  accepted 
criteria  for  transition  from  SMR  to  CMR  and  CMR  to  DMR.  A  worthwhile 
additional  necessary  criterion  for  the  SMR  to  CMR  transition  was  later  made  by 
Shirouzu  and  Glass  [16]. 

The  necessary  conditions  for  the  transition  from  RR  to  MR  in 
nonstationary  flow  are  still  not  fully  understood.  Work  by  Hornung  and 
Kychakoff  [17]  and  Hornung,  Oertel ,  and  Sandeman  [18]  brought  forth  another 
hypothesis,  the  'sonic'  criterion.  It  is  difficult  to  resolve  experimentally 
from  the  'detachment'  criterion,  and  thus  the  'von  Neumann  paradox'  was  still 
unsolved.  It  was  suggested  by  Hornung  et  al  [18]  and  later  investigated  by 
Hornung  and  Taylor  [19]  that  the  ‘von  Neumann  paradox'  was  due  to  the  viscous 
boundary  layer  which  develops  behind  the  shock  on  the  wedge  surface.  Shirouzu 
and  Glass  [16]  supported  this  theory  but  the  evidence  was  still  inconclusive. 

Hornung  and  Taylor  [19]  investigated  the  influence  of  viscosity  on  the 
RR-^-^MR  transition  boundary  by  doing  experiments  in  argon.  They  held  Mach 
number  fixed,  and  varied  the  wedge  angle  to  see  the  behavior  of  the  size  of 
Mach  stem  in  MR.  This  was  done  for  various  initial  pressures.  They  assumed 
transition  to  RR  when  the  Mach  stem  size  approached  zero.  Initial  pressure 
was  found  to  influence  the  position  of  the  transition  boundary. 

It  was  the  intent  of  this  work  to  further  explore  the  effect  of  the 
boundary  layer  on  the  RR*-*>MR  transition  boundary,  although  from  a  different 
approach  than  Hornung  and  Taylor  [19],  The  approach  taken  here  was  to  do 
experiments  in  air  over  a  wide  range  of  Mach  numbers  and  initial  pressures, 
and  to  find  transition  boundary  points  by  doing  both  MR  and  RR  experiments,  in 
close  proximity  to  the  boundary. 


In  this  investigation,  the  objectives  were: 


1.  To  determine  precisely  several  RR*-*»MR  transition  boundary  points 
by  doing  experiments  very  close  to  the  boundary. 

2.  To  use  the  lowest  pressure  possible  to  maximize  viscous  effects. 

3.  To  obtain  the  deviation  from  inviscid  theory  of  the  reflected  wave 
angle  in  RR  at  the  RR**-*TiR  transition  boundary. 

4.  To  predict  the  RR-**>MR  transition  boundary  in  air  for  any  initial 
pressure  and  incident-shock-wave  Mach  number. 

The  present  research  continues  a  long  program  of  study  of 
oblique-shock-wave  reflection  supervised  by  Professor  I.  I.  Glass  at  the 
University  of  Toronto  Institute  for  Aerospace  Studies  (UTIAS).  Previous 
researchers  include  Molder  [20],  Weynants  [21],  Law  and  Glass  [14], 
Ben-Dor  [22],  Ben-Dor  and  Glass  [15,23],  Ando  [24],  Ando  and  Glass  [25],  Lee 
and  Glass  [26],  Shirouzu  and  Glass  [16],  Hu  and  Shirouzu  [27],  Hu  [28], 
Deschambault  [29],  Deschambault  and  Glass  [30],  Hu  and  Glass  [31,32]  and  Li 


Chapter  2 

DESCRIPTION  OF  OBLIQUE-SHOCK-WAVE  REFLECTION 


In  this  chapter,  the  types  of  reflection  and  conditions  associated  with 
each  are  discussed  for  inviscid  flow.  The  effects  of  viscosity  on  shock-wave 
reflection  are  treated  in  Chapter  3. 


2.1  Dynamics  of  Oblique-Shock-Wave  Reflection 


When  a  shock  wave  passes  over  a  wedge  in  a  shock  tube,  two  simultaneous 
processes  occur  [15]: 

(a)  deflection  of  the  flow  at  the  wedge  corner,  and 

(b)  reflection  of  the  incident  shock  wave  from  the  wedge  surface. 

The  first  of  these,  deflection  at  the  corner,  can  take  place  via 

supersonic  turning  through  an  attached  oblique  shock  wave  (Fig.  6a, b),  or 
subsonic  turning  with  a  detached  shock  wave  (Fig.  6c).  In  the  attached  shock 
wave  cases,  the  shock  will  be  straight  if  the  flow  behind  it  is  supersonic 
(Fig.  6a),  or  curved  if  the  flow  behind  it  is  subsonic  (Fig.  6b).  The 

necessary  conditions  for  each  configuration  are  obtained  from  gasdynamics[34] . 

The  reflection  process  can  take  place  via  RR,  SMR,  CMR,  or  DMR.  The 
speed  of  the  incident  shock  wave  is  constant.  While  there  is  some  attenuation 
of  the  shock  wave  with  time  due  to  the  wall  boundary  layer  [36],  it  is  not 
significant  over  very  short  wedge  distances.  The  flow  is  pseudostationary, 
and  consequently  a  frame  of  reference  can  be  attached  to  the  reflection  point 

P  in  RR  and  the  triple  point  T  in  MR.  It  can  also  be  made  self-similar  by 

transforming  the  independent  variables  from  x,y,  and  t  to  x/t  and  y/t  [35]. 
Self-similarity  makes  the  shock  pattern  and  the  flow-property  distributions 
retain  identical  shapes,  but  increase  in  size  linearly  with  time. 


•  >■  .vr 


2.2.1  Regular  Reflection  (RR) 

Regular  reflection  is  shown  schematically  in  Fig.  7.  It  is  the  simplest 
of  the  four  types  of  reflection  since  it  has  only  an  incident  shock  wave  I  and 
a  reflected  shock  wave  R.  The  flow  at  reflection  point  P  must  remain  parallel 
to  the  wedge.  Therefore,  the  flow  deflection  through  R  must  be  equal  and 
opposite  to  the  deflection  through  I.  This  can  be  expressed  as: 

6j  -  02  =  0  (2.1) 

There  are  actually  two  possible  solutions  for  the  reflected-shock-wave  angle, 
a  strong  and  weak  solution  [34].  The  weak  solution  is  the  only  one  which  is 
usually  found  in  practice.  The  simplicity  of  the  shock  configuration  makes 
calculation  of  the  flow  properties  fairly  easy.  Appendix  A  contains  a 
description  of  the  method  of  solution  used  by  Hu  and  Shirouzu  [27]. 

The  reflected  wave  can  be  straight  or  curved  near  P  depending  upon  the 
Mach  number  of  the  flow  in  region  2.  If  the  flow  in  region  2  is 
supersonic,  M^  >  1,  relative  to  P,  the  reflected  wave  must  be  straight  since 
disturbances  from  the  wedge  corner  can  not  propagate  toward  P  (Fig.  8a).  For 
subsonic  flow,  M^  <  1,  in  region  2  relative  to  P,  R  will  be  curved  (Fig.  8b) 
since  disturbances  are  propagated  from  the  corner  to  affect  it. 

2.2.2  Mach  Reflection  (MR) 

Mach  reflection  is  characterized  by  the  appearance  of  a  Mach  stem  M,  and 
slipstream  S  (see  Figs.  9  through  11).  The  necessary  condition  that  flow 
remain  parallel  to  the  wedge  surface  implies  that  the  Mach  stem  must  be 
perpendicular  to  the  wedge  surface.  The  Mach  stem  often  has  some  curvature 
[16],  but  is  usually  approximated  as  being  straight  for  some  analyses. 


Consider  the  flow  in  the  vicinity  of  triple  point  T.  The  flow  can  follow 
two  paths;  the  first  through  regions  0,  1,  and  2,  and  the  second  through 

regions  0  and  3.  Regions  2  and  3  are  separated  by  a  contact  surface  or 

slipstream  S.  Across  S,  the  pressures  and  flow  direction  must  be  equal,  and 

the  density,  temperature,  and  entropy  can  differ.  The  boundary  conditions  are 
therefore: 

0j  -  e2  =  e3  (2.2) 

P2  =  P3  (2.3) 

As  mentioned  previously,  the  flow  at  T  can  be  considered 

pseudostationary,  thus  x  remains  constant  as  the  shock  travels  up  the  wedge. 
The  solution  of  MR  is  somewhat  more  difficult  than  for  RR.  Appendix  A 
outlines  the  method  of  solution  by  Hu  and  Shirouzu  [27],  which  also  contains 
RR  and  MR  solutions  for  a  wide  range  of  Mach  number  and  wedge  angle  for 
various  gases. 

The  unique  properties  of  each  type  of  Mach  reflection  are  now  described. 

2.2.2. 1  Single-Mach  Reflection  (SMR) 

A  schematic  diagram  of  SMR  is  shown  in  Fig.  9.  In  SMR,  unlike  CMR  and 
DMR,  the  reflected  shock  R  is  continuously  curved  from  the  triple  point  T  to 
the  wedge  corner  since  the  flow  in  region  2  is  subsonic  relative  to  T. 

2. 2. 2. 2  Complex-Mach  Reflection  (CMR) 

Figure  10  shows  a  schematic  diagram  of  CMR,  which  differs  from  SMR  due  to 

the  appearance  of  a  kink  K  in  the  reflected  shock.  The  kink  arises  from  a 

band  of  compression  waves  which  converge  on  to  the  reflected  shock  at  K  [12]. 
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The  flow  in  region  2  is  supersonic  relative  to  T  for  CMR  and  DMR,  and  the 
reflected  shock  there  is  straight. 


2. 2. 2. 3  Double-Mach  Reflection  (DMR) 

Figure  11  shows  a  schematic  diagram  of  DMR,  which  differs  from  CMR  due  to 

the  appearance  of  an  additional  Mach  shock  (stem)  M'  and  slipstream  S'  at  a 

new  triple  point  T1  (what  used  to  be  the  kink  K  in  CMR).  The  pattern  of  shock 
waves  and  slipstream  is  the  same  for  both  T  and  T‘.  By  applying  the 
gasdynamics  equations  in  the  same  manner  as  for  the  solution  at  T,  the  flow 

velocities  and  states  can  also  be  solved  for  at  T1  [27]. 

2.3  Transition  Criteria 

2.3.1  RR-***MR  Transition 

Several  different  criteria  have  been  proposed  for  the  transition  from  RR 
to  MR.  Each  of  these  is  now  described, 

2.3. 1.1  'Detachment'  Criterion 

The  'detachment'  criterion  was  proposed  by  von  Neumann  [5,6].  It  stated 
that  transition  would  occur  when  the  flow  deflection  through  the  incident 
shock  exceeded  the  maximum  deflection  possible,  0^,  through  the  reflected 
shock  (Fig.  7). 


9i  -  %  ■  0  (2-4) 

At  this  point,  the  boundary  condition  of  tangential  flow  at  the  wedge  surface 
could  no  longer  be  satisfied  with  a  RR  or  'two-shock'  configuration,  and  MR, 
a  'three-shock'  configuration,  would  be  necessary.  This  boundary  is  shown  in 
Fig.  12. 
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Consider  the  shock  polar  [34]  in  Fig.  13.  The  shock  polar  is  for  a 
constant  Mach  number  Mq  in  region  0.  The  'incident  shock1  curve  represents 
all  possible  flow-pressure  ratios  in  region  1  arising  from  a  flow  deflection 
o j  through  an  incident  shock.  Curves  R1  to  R5  represent  possible 
flow-pressure  ratios  in  region  2  after  a  deflection  through  a  reflected 
shock  R.  The  'Mach  stem'  curve  represents  possible  states  in  region  3  arising 
from  the  flow  through  a  Mach  stem  in  MR.  If  any  of  the  curves,  R1  to  R5, 
intersect  the  ordinate  axis,  RR  is  possible  since  the  condition  of  no  net 
deflection  at  the  wedge  surface  is  satisfied  (Eq.  2.1).  If  any  of  the  curves, 
R1  to  R5,  intersects  the  'Mach  stem'  curve,  then  MR  is  possible,  since  the 
necessary  conditions  of  equal  pressure  and  flow  direction  in  regions  2  and  3 
are  satisfied  (Eqs.  2.2  and  2.3). 

Five  separate  cases  are  represented  by  R1  to  R5: 

Rl:  Only  RR  is  possible  --  point  a 
R2:  RR  and  MR  possible  --  point  b 

No  net  deflection,  therefore  X  is  equal  to  zero 
for  MR  (Fig.  9) 

R3:  RR  (pt.  f)  and  MR  (pt.  c)  both  possible 
R4:  RR  (pt.  g)  and  MR  (pt.  d)  both  possible 
Limit  of  RR  reached 
R5:  Only  MR  is  possible  --  point  e 

The  'detachment'  criterion  states  that  as  the  wedge  angle  is  decreased, 
the  flow  will  follow  path  a-b-f-g  (RR)  and  then  jump  back  to  d  (MR).  In 
effect  it  states  that  as  wedge  angle  is  decreased,  with  Mq  constant,  RR  will 
be  maintained  until  it  is  no  longer  physically  possible,  and  only  then  will  MR 
occur.  It  should  be  noted  that  the  jump  from  g  to  d  has  a  sudden  drop  in 
pressure. 
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2.3. 1.2  'Mechanical  Equilibrium1  Criterion 


The  'mechanical  equilibrium1  criterion,  which  was  proposed  by  Henderson 
and  Lozzi  [12],  is  shown  in  Fig.  12  for  frozen  air,  They  investigated  the 
RR-**»MR  transition  boundary  in  unsteady  flow  by  reflecting  shock  waves  from 
concave  and  convex  wedge  surfaces.  During  transition,  both  from  RR  to  MR  and 
MR  to  RR,  neither  unsteady  expansion  waves  nor  compression  waves  were  seen. 
Some  hysteresis  was  noted  depending  on  whether  transition  was  from  RR  to  MR  or 
MR  to  RR.  Henderson  and  Lozzi  [12]  argued  that  the  absence  of  unsteady 
expansion  or  compression  waves  implied  that  the  transition  process  was  steady, 
in  'mechanical  equilibrium'.  This  criterion  for  transition  was  extended  to 
the  pseudosteady  case,  even  though  their  experiments  with  flat  wedge  surfaces 
showed  persistence  of  RR  below  the  'detachment'  criterion  boundary.  These 
points  which  persisted  into  the  MR  region  were  presumed  to  be  MR  with  Mach 
stems  so  small  that  they  could  not  be  seen. 

On  the  shock  polar  diagram  (Fig.  13)  the  'mechanical  equilibrium' 
criterion  is  represented  by  path  a-b-c-d-e.  Beginning  at  point  a,  as  wedge 
angle  is  decreased  the  solution  will  continue  to  be  one  of  RR.  When  point  b 
is  reached  however,  both  RR  and  MR  are  possible.  According  to  the  criterion, 
points  lying  immediately  on  either  side  of  the  transition  boundary  should  have 
the  same  f low-pressure  ratio.  Therefore,  as  wedge  angle  is  further  decreased, 
MR  occurs  and  the  flow  is  represented  by  points  between  b  and  c.  As  wedge 
angle  is  further  decreased,  the  path  goes  to  d  and  e. 

At  transition  (point  b): 


°1  "  °2  °3  0 


(2.5) 


By  comparison,  the  'detachment'  criterion  states  that  at  b,  when  the  wedge 
angle  is  further  reduced,  the  flow  remains  as  RR  until  point  g  is  reached. 
Further  reduction  in  wedge  angle  causes  'detachment'  to  occur  and  there  is  an 
immediate  jump  to  point  d,  along  with  a  sudden  decrease  in  pressure  ratio. 

Many  experimenters  have  shown  RR  to  persist  below  the  'detachment* 
criterion  boundary,  which  is  well  away  from  the  'mechanical  equilibrium' 
boundary  (Fig.  12).  Some  of  Smith's  data  [8],  which  are  quite  extensive, 
clearly  show  the  persistence  of  RR  (Fig.  14).  Some  more  recent  data  from 
Deschambault  [29],  which  cover  a  wider  range  of  Mach  number,  are  shown  in 
Fig.  15.  It  is  the  'detachment'  criterion  that  is  now  universally  accepted 
for  pseudostationary  oblique-shock-wave  reflection. 

2.3. 1.3  'Sonic'  Criterion 

Hornung  and  Kychakoff  [17]  and  Hornung,  Oertel ,  and  Sandeman  [18] 
proposed  yet  another  theory,  the  'sonic'  criterion.  From  a  physical  point  of 
view,  they  suggested  that  for  a  Mach  stem  to  form,  a  length  scale  would  have 
to  arise  as  a  result  of  signals  transmitted  from  the  wedge  corner.  For  this 
to  occur,  the  flow  in  region  2  relative  to  P  in  RR  must  be  subsonic.  This 
transition  criterion  is  best  expressed  as: 

Oj  -  (>2S  =  0  (2.6) 

The  boundary  for  this  criterion  is  very  close  to  that  predicted  by  the 
'detachment'  criterion  since  the  flow  detachment  angle  o  and  the  sonic  angle 
differ  by  only  one  or  two  degrees  [34],  and  a  resolution  between  the  two 
is  difficult  to  achieve  experimentally. 

Hornung  et  al  [18]  also  hypothesized  that  the  persistence  of  RR  beyond 
the  'sonic'  criterion  boundary  was  due  to  the  viscous  boundary  layer  that  is 
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induced  behind  the  shock  waves.  Hornung  and  Taylor  [19]  later  demonstrated 
that  the  transition  boundary  was  affected  by  changes  in  initial  pressure,  in 
a  relationship  consistent  with  boundary-layer  theory.  However,  their 
experiments  were  only  for  one  particular  Mach  number  in  argon.  A  greater 
discussion  of  their  work  can  be  found  in  Chapter  3. 

For  the  case  of  frozen  air,  a  computer  program  was  written  to  find  the 
RR-^-^MR  transition  boundary  for  any  of  the  three  transition  criteria.  This 
program  is  listed  in  Appendix  B.  The  program  contains  provision  for  viscous 
effects,  and  can  be  easily  modified  to  solve  for  any  frozen  gas.  Its  small 
size  and  simple  construction  makes  it  adaptable  to  microcomputers. 

2.3.2  SMR-^-KIMR  Transition 

The  existence  of  CMR  was  first  noticed  by  White  [9]  to  occur  when  the 
Mach  number  in  region  2  relative  to  the  triple  point  was  greater  than  one. 
This  implies  that  transition  may  occur  when: 

Mo  =  1  (2.7) 
£T 

A  mechanism  for  the  transition  was  proposed  by  Henderson  and  Lozzi  [12]  in 
terms  of  a  band  of  compression  waves.  Law  and  Glass  [14]  and  Ben-Dor  and 
Glass  [15]  verified  that  transition  occurs  according  to  Eq.  2.7.  Law  and 
Glass  [14]  also  predicted  the  location  of  the  kink  with  the  equation: 

1~T  /  =  pi  /  Pq  (2.8) 

where  L^  /  L^  is  the  ratio  of  the  horizontal  distance  from  the  wedge  corner  of 
T  and  K  respectively  (Fig.  16).  Bazhenova  et  al  [13],  Ando  and  Glass  [25], 
and  Deschambault  [29]  found  that  this  expression  was  not  valid  near  the 
RR^-W4R  boundary.  They  showed  experimentally  that  the  ratio  approached  unity. 


An  additional  necessary  criterion  was  added  by  Shirouzu  and  Glass  [16]; 
that  the  angle  between  the  incident  and  reflected  shocks,  <s,  must  exceed  90 


degrees. 

6  >  90°  (2.9) 

This  has  been  verified  experimentally  by  Hu  [28]. 

2.3.3  CMR^-*-0MR  Transition 


Henderson  and  Lozzi  [12]  hypothesized  that  transition  from  CMR  to  DMR 
occurred  when  compression  waves  converged  to  form  the  second  Mach  stem,  M1 
(Fig.  11).  Supporting  this  idea.  Law  and  Glass  [14]  and  Ben-Dor  and  Glass 
[15]  showed  experimentally  that  transition  occurred  when  the  Mach  number  in 
region  2  relative  to  the  kink  equalled  one: 


M 


2 


K 


=  1 


(2.10) 
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Chapter  3 

influence:  of  vibrational  excitation  and  viscosity 

3 . 1  Vibrational  Excitation 

At  room  temperatures,  air  behaves  as  a  frozen  gas.  Only  translational 
and  rotational  degrees  of  freedom  are  excited.  For  this  frozen  or  perfect 
state,  the  specific  heat  ratio,  y,  is  constant,  and  the  Rankine-Hugoniot 
relations  (Ref.  34)  can  be  applied  for  flow  across  a  shock  wave. 

At  high  temperatures,  occurring  behind  a  strong  shock  wave,  real-gas 
effects  come  into  play.  Vibrational  excitation,  dissociation,  and  ionization 
may  be  present.  Dissociation  and  ionization  occur  only  for  very  high  Mach 
numbers  [26],  well  outside  the  range  of  experiments  in  this  study.  As  such, 
they  will  not  be  considered  further. 

To  determine  the  effect  of  vibrational  excitation,  two  things  must  be 
considered;  the  magnitude  of  excitation  at  equilibrium,  and  the  rate  at  which 
equilibrium  is  approached. 

3.1.1  Effect  on  the  RR-**-MR  Boundary 

At  equilibrium,  the  extent  of  vibrational  excitation  depends  only  on 
temperature  (see  App.  C).  However,  the  relaxation  length  to  equilibrium 
depends  on  the  temperature  and  pressure  of  the  gas.  Oxygen  is  excited  at  a 
lower  temperature  than  nitrogen.  The  effect  of  excitation  of  either 
constituent  on  the  RR-<-^MR  boundary  is  to  lower  it,  as  shown  in  Fig.  17.  It 
can  be  seen  that  the  maximum  effect  of  oxygen  excitation  occurs  at  M^  ~  6,  and 
the  maximum  effect  of  oxygen  and  nitrogen  excitation  occurs  at  ~  10.  It 
was  approximated  that  air  consists  of  20  oxygen  and  80%  nitrogen.  A  more 
detailed  description  of  vibrational  relaxation  is  contained  in  Appendix  C. 


For  vibrational  equilibrium,  the  flow  across  a  shock  wave  does  not  have 
a  closed-form  solution  as  does  a  frozen  flow.  It  must  be  solved  numerically 
using  the  equations  of  continuity,  momentum,  energy,  and  the  equation  of 
state. 

3.1.2  Role  of  Relaxation  Length 

Upon  passing  through  a  shock  wave,  the  translational  and  rotational 
degrees  of  freedom  are  fully  excited  within  a  few  mean  free  paths,  the 
thickness  of  the  shock.  The  vibrational  degrees  of  freedom  take  much  longer 
to  reach  equilibrium.  The  size  of  the  relaxation  zone  is  dependent  on 
temperature  and  pressure  (App.  C). 

To  determine  whether  an  assumption  of  frozen  or  equilibrium  flow  is 
appropriate,  the  relaxation  length  must  be  compared  to  a  characteristic  length 
of  the  region  being  studied  behind  the  shock  wave.  If  the  relaxation  length 
is  much  greater  than  the  characteristic  length,  then  the  flow  in  the  region  of 
interest  is  in  essence  frozen.  Conversely,  if  the  relaxation  length  is  much 
smaller  than  the  characteristic  length,  then  an  assumption  of  equilibrium  flow 
is  more  appropriate.  Flows  are  considered  in  nonequilibrium  where  the 
relaxation  length  is  of  comparable  size  to  the  characteristic  length,  and  must 
be  treated  in  a  more  precise  way. 

3.2  Influence  of  Viscosity 

3.2.1  Characteristics  of  Shock-Wave-Induced  Boundary  Layer 

When  a  shock  wave  passes  over  a  surface,  it  induces  a  velocity  in  the  air 
behind  the  shock  [36].  A  boundary  layer  grows  behind  the  shock  wave,  due  to 
the  friction  between  the  moving  air  and  the  stationary  wedge  surface.  This 
boundary  layer  has  a  significant  effect  on  the  flow  at  the 


shock-wave-reflection  point  P.  In  a  laboratory-f ixed  frame  of  reference,  the 

flow  velocity  at  the  wall  must  be  zero  (Fig.  18a).  It  is  convenient  though, 

to  attach  the  reference  frame  to  the  point  where  the  shock  intersects  the 

surface  (Fig.  18b).  In  this  reference  frame,  the  flow  velocity  at  the  wall 

must  be  that  of  the  wall,  relative  to  the  shock.  It  should  be  noted  that  in 

this  new  shock-fixed  reference  frame,  the  boundary-layer  displacement 
★ 

thickness,  .5  ,  has  a  negative  value,  or  the  wall  acts  as  a  sink.  Instead  of 
impeding  the  flow,  the  shear  stress  at  the  wall  tends  to  'sweep'  the  flow  away 
from  the  shock  wave. 

The  boundary-layer  profile  can  be  calculated  by  solving  the  transformed 
laminar-boundary-layer  equations  outlined  by  Mirels  [37,38],  It  will  be  shown 
later  that  the  assumption  of  laminar  flow  is  appropriate,  based  on  the 
Reynolds  number  of  the  experiments  performed.  A  summary  of  the  equations  used 
and  method  of  solution  is  contained  in  Appendices  D  and  E.  Normalized 
velocity  and  temperature  profiles  are  shown  in  Figs.  19  and  20  for  the 
boundary  layer  behind  a  shock  wave  travelling  at  =  2.0,  at  pQ  =  2  kPa  and 
Tq  =  300  K.  While  the  boundary-layer  size  is  significantly  affected  by 
pressure,  and  to  a  lesser  extent  by  temperature  (see  App.  D),  the  shape  of  the 
profiles  is  scarcely  influenced  by  either.  Very  rapid  changes  in  velocity  and 
temperature  occur  close  to  the  wall.  Half  of  the  total  change  in  velocity  and 
temperature  within  the  boundary  layer  occurs  in  the  first  12%  nearest  the 
wall.  It  is  interesting  to  note  that  the  temperature  and  velocity  profiles 
are  quite  similar,  but  at  the  outer  edge  of  the  boundary  layer,  the  normalized 
temperature  only  reaches  0.98  compared  to  the  normalized  velocity  of  0.99. 

The  most  important  thing  about  the  equations  is  the  dependence  of  the 
various  size  parameters  on  density  or  initial  pressure,  and  the  variation  of 
boundary-layer  size  with  distance  from  the  shock  wave. 
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The  relationships  are: 


*  4 

6,  6  *  x? 


x  x*  -£  ~h 

6,  6  «  p  *  or  p 


(3.1a) 


(3.1b) 


This  implies  that  if  shock-wave  reflection  is  influenced  by  viscous  effects 

for  a  given  gas,  Mach  number  M  and  wedge  angle  0  are  not  sufficient  to 

s  w 

define  the  flow;  initial  pressure  or  density  must  also  be  specified.  While 
viscosity  may  not  have  a  great  influence  where  the  Reynolds  number  is  high 
(high  initial  pressures),  many  shock-tube  experiments  are  done  at  very  low 
initial  pressures  where  viscous  effects  are  important. 


3.2.2  Physical  Model  of  the  Reflection  Process  in  RR  with  a 
Shock-Wave-Induced  Boundary  Layer 


As  noted  before,  many  researchers  have  found  significant  differences 
between  experiment  and  the  'detachment'  criterion  boundary,  the  'von  Neumann 
paradox'  (Figs.  14  and  15).  While  this  difference  is  one  indication  of 
possible  viscous  effects,  additional  information  is  gained  by  looking  at  the 
reflected  wave  angle  in  RR.  If  the  wave  angles  near  the  reflection  point  P 
are  measured  off  photographs  from  RR  experiments,  and  the  flow  deflections 
through  them  calculated  using  the  Rankine-Hugoniot  equations,  it  is  found  that 
the  flow  deflection  through  the  reflected  wave  R  is  not  equal  and  opposite  to 
that  through  the  incident  wave  I.  This  difference  in  angles  is  defined  as  the 
displacement  wedge  angle,  and  is  written  as: 


°1  -  02  =  °d 


(3.2) 


and  is  illustrated  in  Fig.  21.  Displacement  wedge  angle,  0^,  is  a  measure  of 
the  deviation  of  RR  experiments  from  inviscid  theory. 


Reflected  wave  angles  u>'  from  RR  experiments  by  Ben-Dor  [22]  and 

Deschambaul t  [29]  were  obtained  by  remeasuring  the  angle  6  between  the 

incident  and  reflected  shock  waves  at  a  distance  of  1  mm  from  the  reflection 

point  P,  and  the  results  for  displacement  wedge  angle  are  shown  in 

Figs.  22  and  23  for  ew  =  50s  and  60°  respectively.  Error  in  measurement 

should  be  quite  small  since  the  reflected  wave  R  is  straight  in  the  vicinity 

of  the  reflection  point  P  for  these  cases.  Figure  22,  which  is  for  RR,  has 

several  theoretical  curves  corresponding  to  6^  =  0°,  1°,  and  2  .  The  curve 

for  0^  =  0  extends  only  up  to  M$  =  1.8,  since  for  Mach  numbers  greater  than 

this  value  and  @w  =  50°,  RR  solutions  do  not  exist.  Similarily,  for  =  T , 

there  are  no  RR  solutions  for  e  =  50°  in  the  range  2.1  <  M  <  3.4.  In 

w  s 

Fig.  23,  the  curves  cover  the  whole  range  of  Mach  number  since  RR  solutions 

exist  for  all  M  when  e  =  60°. 

s  w 

In  Fig.  22,  where  there  is  a  great  deal  of  data,  the  displacement  wedge 
angle  ed  can  be  seen  to  increase  as  Mach  number  increases.  In  the  most 

extreme  case,  for  M$  =  7.3  and  @w  =  50°,  the  displacement  wedge  angle  ed 
exceeds  3C.  Concurrent  with  an  increase  in  Ms,  the  experiments  had  lower 
initial  pressures,  and  it  is  believed  that  this  was  the  cause  of  the  higher 
deviations.  It  is  expected  that  the  same  phenomenon  will  influence  the 

transition  boundary.  Figure  23  also  shows  a  deviation  from  the  inviscid 

prediction,  but  insufficient  data  makes  trends  indeterminate. 

In  a  general  sense,  for  both  RR  and  MR,  the  boundary  layer  serves  as  a 
sink.  Flow  is  drawn  away  from  the  region  by  the  negative  displacement 
thickness.  If  the  negative  thickness  is  added  to  the  wedge  surface  to  give  a 
displaced  wall,  and  the  result  treated  as  an  inviscid  approximation  to  the 
viscous  flow  field,  a  schematic  representation  of  the  physical  processes 
involved  is  obtained  [19].  This  is  illustrated  in  Fig.  24.  According  to  the 


model,  an  expansion  fan  exists  at  the  foot  of,  and  behind  the  incident  shock 

wave,  and  compression  waves  subsequently  coalesce  to  form  the  reflected  shock 

wave.  The  compression  waves  are  due  to  the  decrease  in  rate  of  growth  of  the 

★ 

displacement  thickness  <5  .  The  angle  at  which  flow  enters  the  boundary  layer 
is  henceforth  referred  to  as  'boundary-layer  flow  entry  angle',  6^ ,  and  its 
distribution  with  respect  to  x  is  of  the  form  (see  App.  D): 

tan(eb-|)  =  1  (6*/x^)  x"^  =  x-^  (3.3) 

2 

where  is  a  constant  based  on  the  flow  properties  in  the  freestream, 
external  to  the  boundary  layer.  The  variation  of  with  axial  distance,  as 
defined  by  Eq.  3.3,  is  shown  in  Fig.  25  for  a  typical  boundary  layer  in  air. 
At  x=0,  e..  =  90°  and  the  flow  enters  the  boundary  layer  normal  to  the  surface 
(see  Fig.  24).  At  x  =  1mm,  6^  is  nearly  0°,  and  the  flow  is  nearly  parallel 
to  the  surface.  While  this  model  gives  a  qualitative  illustration  of  the  flow 
processes  involved,  quantitative  information  in  the  immediate  vicinity  of  the 
reflection  point  P  can  not  be  obtained  from  the  simple  model  for  the  following 
reason.  The  expansion  fan  right  at  point  P,  which  turns  the  flow  normal  to 
the  wedge  surface,  is  not  an  accurate  representation  of  the  real  process,  but 
the  result  of  a  singularity  in  the  model  which  occurs  at  x  =  0  (see  Eq.  3.3) 
Since  the  distances  being  dealt  with  are  so  small,  one  must  consider  the  mean 
free  path  of  the  gas  molecules,  and  whether  or  not  the  flow  still  behaves  as 
a  continuum.  As  an  example,  for  a  RR  at  M$  =  2.0,  Pp  =  2.0  kPa,  and 

_5 

Tp  =  300  K,  the  mean  free  path  of  air  is  calculated  to  be  2.5  x  10  mm  based 
on  region  1  flow  Since  the  concept  of  a  continuum  (which  the  model  assumes) 
can  only  be  applied  over  a  scale  greater  than  several  mean  free  paths,  it  is 
unwise  to  believe  any  quantitative  results  within  this  distance  of  P.  The 
expansion  fan  and  subsequent  compression  waves  within  this  distance  of  P  are 
therefore  probably  quite  different  than  what  is  drawn  in  the  model.  Accurate 
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modelling  in  this  region  would  require  kinetic  theory  which  is  beyond  the 
scope  of  this  work.  Note  also  that  the  size  of  these  features  is  so  small 
that  they  can  not  be  seen.  The  model  is  not  based  on  what  is  known  to  occur, 
but  what  is  believed  will  occur. 

It  is  quite  clear  from  the  model  that  the  reflection  region  is  very 
complex.  To  calculate  the  growth  of  the  boundary  layer,  the  freestream 
conditions  outside  the  boundary  layer  must  be  known.  These  change  with 
distance  from  the  reflection  point  however,  ranging  from  region  1  flow 
properties  just  behind  the  incident  shock  wave,  to  region  2  flow  properties 
further  downstream.  As  a  simplification,  the  boundary-layer  growth  is 
calculated  assuming  constant  freestream  conditions  based  on  region  1  flow 
properties. 

The  results  for  boundary-layer  size  at  a  distance  of  1  mm  from  the 
reflection  point  P  are  shown  in  Fig.  26.  As  anticipated,  the  displacement 
thickness  <5*  is  zero  for  M$  =  1  since  there  is  no  flow  perturbation  behind  a 
Mach  wave,  and  6*  levels  off  at  high  Ms. 

According  to  the  prediction  for  boundary-layer  growth,  the  boundary-layer 
flow  entry  angle  0^  will  range  from  90°  at  point  P,  to  0  at  x  =  °°. 
However,  only  a  small  region  near  P  should  influence  the  reflection  process. 
The  size  of  this  region  is  the  characteristic  length  of  the  boundary  layer. 
In  the  case  of  RR,  the  characteristic  length  is  taken  to  be  the  distance  from 
the  reflection  point  P  at  which  the  experimentally  determined  displacement 
wedge  angle  6^,  matches  the  predicted  boundary-layer  entry  angle  0^ •  This 
characteri Stic  length  can  be  calculated  by  setting  equal  to  0^  in 

Eq.  3.3,  and  solving  for  x.  This  definition  of  characteristic  length  is 
somewhat  arbitrary.  It  is  not  expected  that  something  unique  occurs  where  0^ 
and  o^i  match,  but  the  distance  x  at  which  this  happens  is  believed  to  be 


representative  of  the  scale  of  the  region  of  the  boundary  layer  which  affects 
the  reflection  process. 

3.2.3  Effect  of  Shock-Wave-Induced  Boundary  Layer 
on  the  RR-**-MR  Transition  Line 

Henderson  and  Lozzi  [12]  performed  shock-wave-reflection  experiments  with 
a  symmetrical  wedge  in  unsteady  flow  to  eliminate  any  sidewall  boundary-layer 
effects,  and  demonstrated  that  either  the  'sonic'  or  'detachment'  criteria  was 
correct.  As  noted  previously,  it  is  very  difficult  to  experimentally  resolve 
the  difference  between  the  two  criteria.  Hornung  and  Taylor  [19]  did  an 
interesting  set  of  experiments  in  argon  to  see  the  effect  of  initial  pressure 
on  the  RR*-*-MR  transition  boundary.  Using  a  variable-angl e  wedge,  they  did 
four  series  of  experiments,  each  with  a  different  initial  pressure.  The  Mach 
number  was  held  constant.  In  each  series,  the  wedge  angle  was  varied,  and 
they  measured  the  size  of  the  Mach  stem  in  MR.  The  point  at  which  the  Mach 
stem  would  vanish  was  extrapolated  from  each  series,  and  taken  to  be  the  point 
where  transition  would  occur.  The  results  of  the  four  series  showed  that  the 
rr-*h^MR  transition  boundary  was  shifted  downward.  The  magnitude  of  the  shift 
was  greater  when  the  initial  pressure  was  lower,  in  a  manner  consistent  with 
Eq.  3.1b.  Unfortunately,  their  results  were  only  for  one  shock-wave  Mach 
number,  M$  =  5.5. 

As  seen  in  Figs.  22  and  23,  the  boundary  layer  alters  the  boundary 
condition  of  no  net  flow  deflection  at  the  wall  (Eq.  3.2).  To  predict  a  new 
transition  boundary,  the  'detachment'  criterion  might  still  be  used,  though  in 
a  modified  form.  Recall  that  this  criterion  stated  in  principle  that 
transition  would  occur  when  RR  was  no  longer  physically  possible.  At  this 
point,  the  maximum  possible  deflection  would  occur  through  the  reflected  shock 
wave.  Using  this  same  criterion,  but  altering  the  boundary  condition  at  the 


wall  (Eq.  2.1),  the  RR^-^MR  transition  boundary  criterion  can  be  written  as: 


The  effect  of  various  values  of  6^  on  the  transition  boundary  is  shown  in 
Fig.  27.  For  each  degree  of  displacement  wedge  angle,  the  boundary  shifts 
downward  by  approximately  0.7e.  The  shift  in  boundary  varies  linearly  with 


6^  over  the  range  shown.  The  distribution  of  o'  along  the  RR-w-MR 
transition  boundary  is  shown  in  Fig.  28  for  various  values  of  e^.  The 
reflected  wave  angle  u>‘  is  not  affected  much  by  e^.  At  high  M  ,  maximum 
shift  of  w'  is  obtained,  a  reduction  of  w'  of  0.44°  for  every  1°  shift  in 
ed.  Like  the  RR-eH^IR  transition  boundary,  the  shift  in  w'  varies  linearly 
with  6^ . 

It  now  remains  necessary  to  experimentally  determine  the  deviation  of  the 
RR-»^MR  transition  boundary  from  the  inviscid  'detachment'  criterion  boundary, 
and  to  correlate  the  results  with  the  predicted  boundary-layer  size.  This  can 
be  done  by: 

(1)  conducting  experiments  to  obtain  several  transition  boundary  points, 

(2)  measuring  the  reflected  wave  angle 

(3)  calculating  the  net  flow  deflection  through  the  incident  and 
reflected  shock  waves  to  get  the  displacement  wedge  angle  o^, 

(4)  equating  this  value  to  the  predicted  distribution  of 
boundary-layer  flow  entry  angle  0^  using  Eq.  3.3  to  get  a 
characteristic  length. 

A  look  at  flow  properties  of  the  inviscid  RR  solution  can  indicate  which 
parameters  should  affect  the  influence  of  the  boundary  layer  on  the  reflection 
process,  or  displacement  wedge  angle  i^.  According  to  the  simplified  model  of 
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Fig.  24,  the  boundary  layer  in  RR  forms  in  region  1,  behind  the  incident  shock 
wave.  In  addition,  the  reflected  shock  wave  is  due  to  the  compressive  effects 
of  the  displaced  wall  which  must  turn  the  flow  in  region  1  back  parallel  to 
this  wall.  The  reflection  process  is  therefore  dependent  on  the  flow 
parameters,  ,  p^,  T^,  in  region  1.  The  boundary  layer  has  been  modelled 
using  the  displacement  thickness  6*,  which  is  a  measure  of  the  increase  of 
mass  flow  in  the  boundary  layer  relative  to  the  freestream.  Displacement 
thickness  is  defined  at  freestream  conditions  M^,  p^,  T^,  and  by  using  it  to 
represent  the  boundary  layer,  the  influence  of  p^  and  T^  are  accounted  for. 
The  effect  of  and  on  the  reflection  process  in  the  presence  of  a 
boundary  layer  is  still  unknown. 

The  variation  of  o,  and  M.  with  M  is  shown  in  Figs.  29  and  30.  It  can 
be  seen  that  the  two  parameters  are  relatively  stationary  at  high  Mach  number 
(M  >  5.0),  and  change  quite  rapidly  at  low  Mach  number  (M  <  3.0).  The 
boundary-layer  displacement  thickness  is  also  fairly  constant  for  high  M$ 
(Fig.  26).  The  three  parameters  which  are  believed  may  affect  the  RR-*-*-MR 
transition  boundary  shift  due  to  viscosity,  M^,  0^,  and  boundary-layer  size 
(for  fixed  p^)  are  all  approximately  constant  at  high  M$.  It  should  therefore 
be  anticipated  that  for  a  fixed  initial  pressure,  the  displacement  wedge  angle 
Qj ,  and  therefore  the  shift  of  the  RR-*-*41R  transition  boundary  due  to 
viscosity,  should  also  be  relatively  invariant  at  high  M  .  As  a  result,  a 
viscous  RR-*^MR  transition  boundary  should  be  quite  flat  at  high  (see 
Fig.  27). 
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Chapter  4 
EXPERIMENTAL  WORK 


4.1  Shock-Tube  Facility 

Experiments  were  carried  out  in  the  UTIAS  10cm  x  18cm  Hypervelocity  Shock 
Tube.  This  facility  has  the  capability  of  both  cold-gas  driven  and 
combustion  driven  runs.  However,  combustion  driven  runs  were  not  attempted 
owing  to  the  possibility  of  damaging  the  test-section  windows  when  operating 
with  high  wedge  angle  [39].  A  full  description  of  the  facility  and  its 
capabilities  can  be  found  in  Boyer  [40]  and  Ben-Dor  and  Whitten  [41]. 

By  increasing  the  pressure  in  the  compression  chamber  on  one  side  of  a 
diaphragm  until  it  ruptured,  shock  waves  were  sent  down  the  channel.  The 
driver  gas  was  added  slowly  so  that  after  the  diaphragm  ruptured,  the  shock 
speed  would  not  be  influenced  during  its  time  of  travel  by  the  additional  mass 
flow.  Shock  speeds  were  controlled  by  varying  the  type  of  driver  gas,  the 
diaphragm  thickness  and  pressure  ratio  across  the  diaphragm.  The  diaphragms 
consisted  of  several  sheets  of  mylar  polyester  stacked  together.  The  sheets 
were  available  in  several  nominal  thicknesses,  and  by  proper  combination 
nearly  any  desired  thickness  could  be  achieved.  There  is  a  practical  limit  to 
the  maximum  thickness  which  can  be  used  since  very  thick  diaphragms  do  not 
break  properly  and  their  breaking  pressure  is  unpredictable.  For  this  reason, 
the  maximum  overall  diaphragm  thickness  used  was  1.07  mm.  The  driver  gases 
used  were  CO^  and  He.  Shock-wave  Mach  numbers  up  to  M£  =  6.5  were  achieved. 
While  H^  as  a  driver  gas  would  give  even  higher  shock  speeds,  the  possible 
risk  of  explosion  made  its  use  undesirable.  For  test  gas  pressure  under  26.7 
kPa  (200  torr),  commercially  available  Medical  Grade  Air  was  used.  All  other 
runs  used  air  from  the  laboratory  as  the  test  gas. 


Initial  gas  temperature  was  measured  by  a  mercury-bulb  thermometer 
embedded  in  the  shock  tube  wall.  A  minimum  of  5  minutes  was  allowed  between 
entry  of  the  test  gas  and  reading  from  the  thermometer  to  allow  the 
temperature  to  stabilize.  The  thermometer  scale  was  marked  at  0.1  C 
intervals. 

The  initial  gas  pressure  was  measured  by  a  series  of  Wallace  and  Tiernan 
type  FA  160  dial  gauges.  The  range  and  maximum  error  of  the  gauges  are  shown 
below. 


Range 


1. 

0 

-  50 

torr 

( 

6.67 

kPa ) 

2. 

0 

-  200 

torr 

( 

26.7 

kPa ) 

3. 

0 

-  400 

torr 

( 

53.3 

kPa) 

4. 

400 

-  800 

torr 

(106.7 

kPa ) 

Absolute  Error 
0.2  torr 
1.0  torr 
2.0  torr 
2.0  torr 


To  guard  against  possible  leakage  of  air  into  the  shock-tube  test  section, 
test  gas  pressure  was  checked  at  the  time  of  admission,  and  also  5  minutes 
after.  Leaks  which  would  cause  a  variation  in  pressure  of  greater  than  0.05 
torr  (6.7  Pa)  at  the  time  of  firing  the  shock  tube  were  unacceptable  and  the 
runs  were  aborted.  This  corresponds  to  a  leak  rate  of  approximately  0.01 
torr/min.  and  was  found  to  be  extremely  rare. 


The  incident-shock-wave  speed  was  determined  by  measuring  the  time 
interval  between  a  common  starting  point  and  several  stations  distributed  down 
the  channel.  Time  intervals  were  measured  by  5  digital  counters  using  a  1  MHz 
oscillator  as  a  common  time  base.  Trigger  signals  from  each  station  were 
obtained  using  Atlantic  LD-25  blast-wave  pressure-transducers .  Start  signal 
for  all  the  counters  came  from  station  D,  and  the  5  counters  were  stopped  by 
signals  from  stations  F,  G,  H,  I,  and  J.  A  schematic  diagram  of  the  set-up  is 
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shown  in  Fig.  31.  To  allow  for  deceleration  of  the  shock  wave,  owing  to  the 
sidewall  boundary-layer  growth,  a  second-order  best-fit  curve  was  applied  to 
the  measurements  from  stations  F,  G,  H,  and  I,  and  extrapolated  to  give  a 
velocity  at  the  center  of  the  test  section.  The  interval  at  station  J  was  not 
used  because  the  shock-wave  speed  changes  significantly  after  it  has  passed 
over  the  wedge. 

Shock-wave-reflection  patterns  were  recorded  using  a  23  cm  dia. 
f ield-of-view  Mach-Zender  i nterferometer .  With  the  exception  of  4  experiments 
where  shadowgraphs  were  taken,  infinite-fringe  interferograms  were  used.  Two 
wavelengths  were  used  for  the  interferograms  (3971  A  and  6943  A)  to  establish 
both  qualitative  and  quantitative  results.  The  3971  A  interferograms  give 
twice  as  many  fringes  as  the  6943  A  interferograms ,  but  suffer  from  lack  of 
contrast.  This  type  of  recording  was  chosen  because  it  gives  a  great  deal  of 
quantitative  density  information,  and  the  presence  of  fringes  makes  it  easier 
to  see  the  Mach  stem,  even  when  it  is  quite  small.  The  4  shadowgraphs  were 
taken  to  compare  resolution  between  them  and  the  interferograms ,  and  were 
found  to  be  inferior. 

The  light  source  for  the  interferometer  was  a  giant-pulse  ruby-laser 
incorporating  a  TRG  model  104A  laser  head,  a  TRG  model  2113-1  harmonic 
generator  and  TRG  Pockels  Q-switch.  The  flashlamp  capacitor  bank  of  the  laser 
and  the  Q-switch  were  fired  at  the  appropriate  times  by  delaying  the  trigger 
signal  from  stations  F  and  I  respectively.  The  capacitor  bank  was  fired  at 
900_s  before  the  Q-switch  in  order  to  get  a  single  giant  pulse  from  the  laser. 
In  this  mode,  the  laser  delivers  0.8  Joules  of  light  energy  in  a  15  ns  pulse. 
Interferograms  were  recorded  on  Kodak  Royal -X  Pan  film,  and  developed  in  Kodak 
HC-110  developer.  The  negatives  were  developed  using  twice  the  recommended 
time,  to  increase  the  effective  film  speed  (nominal  speed  is  1250  ASA  ,  and 


4.2  Range  of  Experiments 


Experiments  were  done  to  determine  several  RR-*  ►MR  transition  boundary 
points.  Equipment  limitations  often  determined  which  initial  pressure  couid 
be  used.  In  the  low  Mach  number  range  (M  <  1.5),  pressures  had  to  be  near 
atmospheric  so  that  the  pressure  jump  across  the  incident  shock  wave  was 
sufficient  to  trigger  the  timing  equipment  and  the  laser.  At  very  high  Mach 
numbers  (M  <  4.0),  very  low  initial  pressures  were  needed  to  achieve 
sufficiently  high  pressure  ratios  across  the  diaphragm,  but  at  the  same  time 
must  be  high  enough  to  trigger  the  equipment.  Low  initial  pressures  are 
extremely  beneficial  since  they  enhance  viscous  effects. 


! 
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Experiments  were  performed  over  a  range  1.1  <  M$  <  6.5  using  wedges  of 
42”,  45c,  47J,  and  48'.  Pressures  ranged  from  400  Pa  (3  torr)  at  high  M  ,  to 

101  kPa  (760  torr)  at  very  low  Typically,  to  determine  a  particular 

transition  boundary  point,  wedge  angle  and  diaphragm  thickness  were  held 
constant.  Runs  were  then  done  at  several  different  initial  pressures.  The 
Mach  number  for  each  run  is  thus  slightly  different  since  the  pressure  ratio 
across  the  diaphragm  is  altered.  A  transition  point  (Pq,  M  ,  o  )  is  thus 
obtained. 

It  is  a  desirable  condition  to  vary  the  initial  pressure,  Pq,  and  shock 
Mach  number,  M  ,  for  a  fixed  wedge  angle,  6^,  particularly  when  looking  for 
a  transition  point  at  high  Mach  number.  As  discussed  in  Chapter  3,  for  a 
constant  pressure  it  is  expected  that  the  RR-<-^MR  transition  boundary  on  the 
(M  ,  w)  plane  is  relatively  flat  at  high  Mg  for  p^  constant. 
Experimentally  then,  if  p^  and  are  fixed,  and  M^  varied,  one  could 

conceivably  be  doing  experiments  parallel  to  the  transition  boundary,  but  not 
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crossing  it.  In  the  low  M  region  where  the  inviscid  transition  boundary  is 
quite  steep,  this  problem  should  not  occur. 

It  was  attempted  in  one  series  of  experiments  (at  low  M  )  to  hold  the 
initial  pressure  constant  and  change  the  diaphragm  thickness  to  control  shock 
speed,  but  it  was  extremely  difficult  to  find  the  proper  combination  of 
diaphragms  to  achieve  the  required  breaking  pressure  and  shock  speed  M  .  This 
method  is  not  recommended.  A  variable-angle  wedge  or  plunger  to  rupture  the 
diaphragm  at  a  desired  pressure  would  solve  most  of  these  problems,  but  the 
UTIAS  shock  tube  is  not  equipped  for  either.  If  further  work  in  finding 
transition  boundary  points  is  to  be  done,  it  is  highly  recommended  that  both 
mechanisms  be  used  to  greatly  simplify  the  experiments. 

4.3  Measurement  of  Data  from  Interferograms 

The  one  drawback  to  recording  shock-wave  reflections  photographically  is 
the  inevitable  problem  of  interpreting  the  data.  Flaws  in  the  optics  and 
grain  of  the  film  limit  the  resolution  of  the  i nterferograms ,  and  can  make 
determination  of  the  type  of  reflection  difficult  and  subjective.  In  the  case 
of  transition  from  RR  to  DMR,  determination  is  rather  simple.  In  DMR  the  two 
Mach  stems  are  separate  and  distinct.  Transition  from  RR  to  SMR  or  CMR  is 
often  difficult  to  see.  The  Mach  stem  in  MR  may  be  so  small  that  its  size  is 
almost  immeasurable,  and  the  fringe  distortion  caused  by  the  slipstream  can  be 
confused  with  that  caused  by  the  boundary  layer. 

One  way  of  separating  the  two  'ls  to  look  at  the  shape  of  the  fringe 
distortion.  A  slipstream  should  slope  toward  the  wedge  surface  and  distortion 
from  it  appear  to  decrease  in  size  as  distance  from  the  Mach  stem  increases. 
A  boundary-layer-induced  distortion,  if  visible,  should  increase  in  size  as 
distance  from  the  reflection  point  increases. 
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Measurement  of  wave  angles  is  also  difficult.  For  optical  reasons,  the 
reflected  shock  wave  appears  to  have  a  finite,  measurable  thickness  (Figs.  39a 
and  39b).  It  also  appears  to  be  quite  'fuzzy'  near  the  reflection  point.  It 
was  found  that  a  precise  and  accurate  measurement  of  angle  could  only  be  made 
at  distances  of  1  mm  or  greater  from  the  reflection  point.  The  angle  between 
the  incident  and  reflected  shock  waves,  6,  was  measured  by  drawing  lines 
tangent  to  the  incident  and  reflected  shock  waves  at  a  distance  of  1mm  from  P, 
and  trigonometrical ly  determining  the  angle.  Both  first  and  second  harmonic 
photographs  were  measured  independently.  If  the  angles  measured  from  the  two 
differed  by  more  than  1.5’,  they  were  remeasured  until  they  were  within  1.5=. 
The  wave  angle  was  taken  as  the  average.  Wave  angles  measured  by  different 
people  using  different  techniques  can  differ  by  as  much  as  3°  and  is  highly 
subjective.  For  this  reason,  one  should  not  be  too  alarmed  by  different  wave 
angles  quoted  for  the  same  pnotograph. 


Chapter  5 

DISCUSSION  OF  RESULTS 


In  this  chapter,  the  experimental  results  will  be  discussed  in  three 
separate  and  distinct  sections: 

(1)  A  qualitative  look  at  some  of  the  features  observed  of  the  4 
types  of  reflection 

(2)  An  analysis  of  the  transition  boundary  shift  obtained,  and  its 
relationship  to  the  boundary  layer 

(3)  A  brief  observation  of  the  effect  of  the  boundary  layer  on 
the  Mach  stem  in  MR 

5 . 1  Observed  Features  RR,  SMR,  CMR  and  DMR 

Figures  39  through  42  illustrate  the  four  types  of  reflection  obtained, 
and  reveal  some  interesting  features  of  the  flow.  Figures  39a  and  39b  show  a 
case  of  regular  reflection  (RR).  No  evidence  of  a  Mach  stem  or  slipstream  is 

seen  at  the  reflection  point  P.  Some  distortion  at  the  wedge  surface 

indicates  the  presence  of  a  boundary  layer. 

There  is  no  physical  evidence  that  the  structure  of  the  reflected  shock 
wave  is  like  the  physical  model  in  Chapter  3  (Fig.  24).  Rapid  curvature  of 
the  reflected  shock  wave  near  the  reflection  point  is  not  seen  (Fig.  39b), 
contrary  to  the  model.  This  would  suggest  that  the  size  of  region  which 

dominates  the  reflection  process  may  be  smaller  than  can  be  seen  in  the 

interferogram,  or  that  the  model  is  inaccurate  and  a  better  one  is  needed. 

For  this  set  of  experimental  conditions,  the  theoretical,  inviscid 
solution  is  that  of  Mach  reflection  (MR),  with  a  triple-point-trajectory  angle 
x  =  2.84  .  This  is  shown  in  Fig.  39a.  The  difference  between  theory  and 
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experiment  is  particularly  notable  in  this  case,  since  the  initial  pressure 


f 

I 


is  quite  high  (p^  =  100.8  kPa).  With  lower  pressure,  one  would  expect  even 
greater  disagreement  between  theory  and  experiment. 

Figures  40a  and  40b  show  a  case  of  SMR.  The  Mach  stem  is  quite  evident, 
as  is  the  slipstream.  Again,  the  Mach  stem  is  much  smaller  than  the  inviscid 
prediction.  Clearly,  the  slipstream  slopes  toward  the  surface,  and  it  is  this 
feature  which  distinguishes  it  from  the  boundary  layer.  In  some  cases  with  an 
extremely  small  Mach  stem,  this  is  the  only  way  to  differentiate  between  RR 
and  MR. 

Figures  41a  and  41b  show  a  case  of  CMR.  The  exact  position  of  the  kink 
K  is  not  certain  due  to  the  apparent  thickness  of  the  shocks,  but  the 
curvature  reversal  is  quite  clear.  It  is  probable  that  some  cases  of  CMR  were 
judged  to  be  SMR.  Near  the  RR^-^MR  boundary,  the  kink  is  very  weak,  and  it  is 
likely  that  in  some  marginal  cases  it  can  not  be  distinguished,  even  though  it 
exists.  The  slipstream  emanating  from  the  triple  point  slopes  down  to 
interact  with  the  boundary  layer  on  the  wedge  surface .  Note  that  this  is  a 
shadowgraph,  so  fringes  do  not  appear.  In  one  aspect,  shadowgraphs  are 
superior.  The  shocks  tend  to  appear  thinner  and  more  well  defined  than  in 
interferograms. 

Finally,  DMR  is  shown  in  Figs.  42a  and  42b.  Again,  as  with  SMR  and  CMR, 
the  triple  point  does  not  follow  the  predicted  trajectory.  In  this  case  the 
difference  is  quite  large  (2.8*),  likely  due  to  the  low  initial  pressure 
(p0  =  1.47  kPa),  and  the  resultant  large  viscous  effects. 

5.2  RR-*-»MR  Transition  Boundary  Shift 

A  summary  of  the  98  experiments  performed  is  contained  in  Table  1.  Tn 
addition,  they  are  plotted  in  the  (M  ,  ■•  )  plane  for  each  wedge  angle  in  Figs. 
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32  to  38.  From  the  series  of  experiments,  several  points  were  selected  which 
were  felt  to  represent  the  limit  of  RR,  transition  boundary  points.  The 
selection  of  transition  points  from  experimental  data  is  often  not  clear  and 
as  a  result  subjective.  The  difficulty  comes  about  in  part  from  the  fact  that 
both  Mach  number  and  initial  pressure  change  in  a  series,  a  limitation 
dictated  by  the  equipment  available  (see  Chapter  4).  Ideally,  for  a  fixed 
wedge  angle,  as  the  initial  pressure  is  dropped  the  Mach  number  will  increase 
gradually.  All  the  experiments  will  be  one  type  of  reflection,  and  then  at 
some  point  the  type  of  reflection  will  change  from  RR  to  MR  or  MR  to  RR.  As 
Mach  number  is  further  increased  the  type  of  reflection  would  not  change  back. 
In  this  ideal  situation  the  point  of  transition  is  quite  distinct.  Remember 
that  a  transition  point  is  defined  by  M  ,  Q^,  and  pg.  Although  far  from 
ideal,  this  behavior  can  be  seen  in  Figs.  32,  33,  35a,  35b,  and  38.  A 
transition  point  is  relatively  easy  to  pick. 

In  Figs.  34,  36,  and  37,  the  transition  point  is  not  so  clear.  The 
regions  of  RR  and  MR  are  not  distinct  from  one  another,  and  there  is  overlap. 
The  progression  of  Mach  number  with  decreasing  initial  pressure  is  not  smooth. 
In  cases  where  MR  experiments  appeared  in  what  was  thought  to  be  a  RR  zone  or 
vice  versa,  runs  were  repeated  in  an  attempt  to  check  their  validity. 
However,  in  these  three  series  of  experiments,  the  RR  and  MR  regions  could  not 
be  clearly  defined.  It  then  became  highly  subjective  in  selecting  a 
transition  point.  In  Fig.  34  there  is  some  overlap  but  it  is  not  too  bad. 
Figure  36  has  quite  a  great  deal  of  it,  and  the  only  portion  where  RR  and  MR 
are  separate  and  distinct  is  near  M$  =  4.0.  Figure  37  has  two  RR  points  in  an 
apparent  MR  zone  (M  =  4.65,  5.28)  but  these  are  not  considered  valid  since 
they  could  not  be  repeated.  Some  of  the  overlap  of  regions  at  low  M^  is 
undoubtedly  due  to  misjudging  of  type  of  reflection.  MR  runs  at  low  M$  are 
very  difficult  to  distinguish  from  RR.  The  points  which  were  judged  to 


represent  the  limit  of  RR,  or  transition,  are  indicated  on  each  of  the  figures 
and  summarized  in  Table  2. 

An  additional  transition  boundary  point  is  taken  from  experiments  by 
Deschambault  [29].  In  a  series  of  experiments  at  (Mg  =  3.71,  0 w  =  49",  Pq  = 
6  kPa)  he  obtained  examples  of  both  RR  and  MR.  It  was  suggested  that  the  RR 
case  was  actually  MR,  its  Mach  stem  not  having  had  time  to  grow  sufficiently 

to  be  seen.  In  the  RR  case  the  shock  had  travelled  60%  of  the  way  up  the 

wedge.  In  the  MR  case,  the  shock  had  travelled  80%  of  the  way,  and  had  a 
triple-point-trajectory  angle  x  =  1.3  .  If  the  RR  case  were  actually  MR  as 
suggested,  then  one  might  expect  that  with  x  =  1.3",  the  Mach  stem  would  be 
approximately  1.4  mm  high  at  60%  up  the  wedge  surface.  A  Mach  stem  of  this 
size  would  be  quite  visible.  A  more  likely  explanation  for  the  difference  is 
that  the  experiment  lies  at  or  very  near  the  RR  MR  transition  boundary. 

From  an  examination  of  Figs.  35a  and  35b,  it  is  quite  clear  that  initial 

pressure  has  an  effect  on  the  position  of  the  transition  boundary.  The  two 

sets  of  runs,  which  have  a  different  range  of  initial  pressures,  have 
significantly  different  Mach  numbers  at  transition.  The  transition  point 
corresponding  to  the  lower  initial  pressure  has  a  higher  Mach  number,  which 
puts  it  further  from  the  inviscid  transition  boundary.  This  is  entirely 
expected  as  viscous  effects  should  be  more  prominent  for  cases  with  lower 
initial  pressure. 

The  analysis  of  the  transition  boundary  points  will  take  the  following 
steps : 

(1)  Comparison  of  the  selected  transition  boundary  points  with 
the  inviscid  'detachment'  criterion  boundary 

(2)  Comparison  of  reflected  wave  angle  w'  with  that  predicted 
by  the  'detachment'  criterion 


(3)  Reduction  of  experimentally  determined  displacement  wedge 


angle  to  a  common  initial  pressure 

(4)  Comparison  of  to  predicted  boundary-layer  size 

(5)  Calculation  of  characteristic  length 

(6)  Prediction  of  a  viscous  RR-*-*»MR  transition  boundary  for  air 

The  selected  transition  boundary  points  are  plotted  against  the  inviscid 

'detachment'  criterion  boundary  in  Figs.  43  and  44.  Notice  that  the  deviation 

from  the  inviscid  boundary  is  mucn  greater  at  high  M$,  due  primarily  to  low 

initial  pressure.  The  highest  deviation  from  the  'detachment'  criterion 

boundary  is  5.4  in  e  ,  which  occurred  for  an  experiment  with  6  =  45 

w  w 

=  6.21,  Pq  =  0.45  kPa.  This  behavior,  persistence  of  RR  well  below  the 
'detachment'  criterion  boundary  at  high  M  ,  is  certainly  not  applicable  to  all 
gases.  Experiments  by  Hu  [28]  in  SFg  and  Ikui  et  al  [42]  in  Freon-12  show 
persistence  of  DMR  above  the  'detachment'  criterion  boundary  at  high 
These  two  gases  are  unique  since  they  both  have  a  low  specific  heat  ratio  y. 
One  should  therefore  be  cautious  in  extending  the  high  M$  results  of  this 
investigation  in  air  to  other  gases  without  doing  a  similar  analysis. 

In  the  low  regime  (Fig.  1-3),  some  of  the  points  seem  quite  close  to 
the  inviscid  boundary,  and  the  difference  might  be  taken  as  insignificant. 
However,  this  is  not  the  case.  Only  the  steep  slope  of  the  boundary  makes  the 
differences  appear  quite  small.  Note  also  (Fig.  26)  that  the  boundary  layer 
should  also  be  much  smaller  at  low  M  . 

It  can  also  be  shown  that  the  'detachment'  criterion,  which  states  that 
transition  will  occur  when  the  maximum  flow  deflection  possible  occurs  through 
the  reflected  shock  wave,  still  seems  to  hold  under  the  influence  of 


are 


viscosity.  In  Fig  45,  measured  values  of  reflected-shock-wave  angle  w' 
compared  to  those  predicted  by  a  modified  'detachment'  criterion  (Eq.  3.4). 
The  variation  of  reflected-shock-wave  angle  u>'  with  displacement  wedge  angle 
6^  as  predicted  by  Eq.  3.4  is  quite  small.  Experimentally  obtained  w' 

values  are  close  to  the  predicted  values,  within  experimental  error  in 
measurement.  This  implies  that  the  concept  of  the  'detachment'  criterion  is 
still  valid. 

Due  to  experimental  limitations,  there  was  a  wide  variety  of  initial 
pressures.  Valid  comparison  of  results  requires  reduction  to  some  standard. 
While  it  may  at  first  seem  desirable  to  reduce  results  to  a  common  Reynolds 
number,  this  is  not  chosen  for  the  following  reasons.  The  boundary-layer  size 
does  not  reduce  to  a  simple  dependence  on  Reynolds  number  (see  App.  D).  Also, 
it  would  be  a  misleading  representation  since  the  only  valid  length  scale, 
which  is  necessary  to  define  a  Reynolds  number,  is  the  characteristic  length 
which  may  not  be  the  same  for  all  the  experiments. 

Instead  it  is  chosen  to  reduce  the  results  to  a  common  initial  pressure. 
From  Eq.  3.1b  it  can  be  seen  that  any  viscous  effect  can  be  scaled  by  initial 
pressure.  Figure  46  displays  the  experimentally  determined  displacement  wedge 
angle  e^,  a  result  of  viscous  effects,  after  reduction  to  an  initial 
pressure  of  2.0  kPa.  Remember  that  is  obtained  by  applying  gasdynamics  of 
frozen  flow  to  calculate  the  flow  deflection  across  the  incident  shock  wave 
and  the  measured  reflected  shock  wave.  As  anticipated,  the  displacement  wedge 
angle  e^,  a  measure  of  departure  from  the  'detachment'  criterion  boundary,  is 
relatively  constant  at  high  M  ,  with  a  value  of  approximately  4°.  At  low  M$, 
0^  is  higher,  near  6.5r.  This  is  an  unexpected  behavior  since  the  predicted 
boundary-layer  size  is  much  smaller  at  low  (Fig.  26).  The  source  of 

scatter  of  the  two  points  near  =  3.8  is  unknown. 
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To  appreciate  this  difference  quantitatively.  Fig.  47  shows  the  ratio 
between  a  predicted  boundary-layer  size  parameter,  and  the  experimentally 


determined  e^.  As  a  measure  of  the  predicted  boundary-layer  size,  the 
boundary-layer  flow  entry  angle  9^  at  a  distance  of  1  mm  from  the  reflection 
point  P  is  used  (see  Eq.  3.3),  and  given  the  special  symbol  n.  This 
boundary-layer  size  parameter,  u  ,  was  chosen  instead  of  6*/  since  it  has 

the  same  sense  as  e^,  both  being  angles.  While  less  fundamental  than  6*  /  y) 
it  is  a  valid  representation,  for  comparitive  purposes,  of  boundary-layer 
size,  since  it  is  proportional  to  &*  /  x^. 


bl 


x  =  lmm 


(5.1) 


The  experimentally  determined  is  much  larger  than  6.  ,  in  the  predicted 
boundary-layer  at  x  =  1  mm.  This  indicates  that  the  characteristic  length  of 
the  boundary  layer  which  dominates  the  reflection  process  is  much  less  than 
1  mm. 


Since  the  displacement  wedge  angle  has  been  determined  experimentally, 
and  the  boundary-layer  size  <5*  /  x^  has  been  predicted,  the  characteristic 
length  can  be  determined  by  putting  the  two  values  into  Eq.  3.3  and  solving 
for  x.  This  has  been  done  and  is  shown  in  Fig.  48.  The  characteristic  length 
ranges  from  0.0025  mm  at  low  to  0.048  mm  at  high  Ms>  With  values  so  small, 
it  would  be  impossible  to  resolve  any  features  in  the  region  of  reflection 
point  P  which  is  dominated  by  viscous  forces.  It  is  also  so  small  that  the 
previous  assumption  of  frozen  flow  and  laminar  boundary  layer  in  the  region  of 
interest  is  correct.  The  smallest  vibrational  relaxation  length  for  the 
experiments  was  found  to  be  33  mm  for  0^  in  Exp.  62  for  which  =  0.45  kPa 
Tq=24.8  C,  Ms  =6.21  (see  Fig.  C-l).  In  addition,  the  maximum  Reynolds  number 

3 

was  of  the  order  of  10  based  on  the  characteristic  length,  much  less  than 
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that  at  which  transition  to  turbulent  flow  (10  )  might  be  expected  to 
occur  [36]. 

The  inability  to  see  any  details  of  the  viscous  reflection  process  near 
P  demonstrates  the  necessity  for  numerical  codes  which  can  simulate  these 
processes.  Inviscid  numerical  simulations  of  shock-wave  reflections  currently 
give  excellent  results  [43],  as  do  codes  which  approximate  the  effects  of 
viscosity  [51],  but  a  solution  of  the  Navier-Stokes  equations  very  close  to 
the  reflection  point  P  is  still  a  step  beyond.  However,  rapid  technological 
advancement  of  supercomputers  should  make  this  possible  in  the  near  future. 

It  was  shown  in  Fig.  47  that  the  ratio  of  experimentally  determined 
displacement  wedge  angle  0^  to  boundary-layer  size  parameter  fi  was  much 
greater  at  low  M  .  In  other  words,  a  boundary-layer  of  a  particular  size  will 
have  much  more  effect  on  RR  at  low  than  at  high  Ms.  It  is  hypothesized 
that  this  ratio  should  be  dependent  on  flow  deflection  through  the 
incident  shock  wave  for  the  following  reason.  The  ability  of  a  boundary  layer 
to  influence  a  reflection  process  should  be  affected  by  the  relative  size  of 
the  forces  involved.  The  effect  of  the  boundary-layer  is  that  of  mass 
transfer  from  the  reflection  point  P.  The  rate  at  which  mass  is  transferred 
toward  the  wedge  surface  after  deflection  through  the  incident  shock  wave  is 
proportional  to  sin(e^).  If  the  flow  deflection  through  the  incident  shock 
wave  is  quite  small,  then  the  influence  of  the  boundary  layer  should  be  quite 
significant  compared  to  the  mass  transfer  towards  the  wall  which  gives  rise  to 
the  reflection  process.  At  high  M  ,  the  mass  transfer  in  the  boundary  layer 
is  less  significant  compared  to  the  mass  transfer  toward  the  wall,  so  it  has 
less  effect  on  RR. 
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Figure  49  shows  the  ratio  <2  /  as  a  function  of  deflection  angle 
through  the  incident  shock  wave.  is  very  close  to  sin(o^).  The 
relationship  appears  to  be  close  to  linear,  as  might  be  expected  for  the 

hypothesis  given  above. 

The  net  result  of  all  this  is  a  prediction  (based  on  experimental 
results)  of  a  viscous  RR-^^MR  transition  boundary  for  air.  If  a  linear 
relationship  between  n  /  and  sin(e^)  is  assumed  (see  Fig.  49),  and  the 
result  applied  to  the  predicted  values  of  and  boundary-layer  size  parameter 

(Eq.  5.1),  a  distribution  of  6^  is  obtained.  When  substituted  into  the 

modified  'detachment'  criterion  (Eq.  3.2)  the  result  is  a  new  RR-*»-MR 

transition  boundary.  Tnis  is  shown  in  Fig. 50,  for  a  nominal  initial  pressure 
of  5.0  kPa.  Scaling  to  any  other  initial  pressure  is  simple,  since  the 
difference  between  the  inviscid  and  viscous  boundaries  is  a  direct  function  of 
initial  pressure  (Eq.  3.1b).  Note  that  Fig.  50  is  quite  similar  to  Fig.  27 
since  the  experimental  values  for  ed  were  found  to  be  constant  at  high  , 
which  is  precisely  what  is  plotted  in  Fig.  27,  curves  of  constant  0^. 

5.3  Behavior  of  X  Near  RR-«»-MR  Boundary 

In  the  series  of  experiments  for  wedge  angles  of  47  and  48  ,  many  MR 
cases  were  obtained,  with  reasonable  variation  of  initial  pressure.  It  is 
expected  that  the  boundary  layer  on  the  wedge  surface  will  act  as  a  sink  and 
'puli'  the  triple  point  down  toward  the  wedge.  The  boundary-layer  size,  and 
presumably  its  influence,  should  increase  with  decreasing  initial  pressure 
(see  Eq.  3.1b).  The  difference  between  the  inviscid  prediction  of 
tri pi e-poi nt-tra jectory  angle  x  and  the  experimentally  obtained  value  should 
therefore  be  affected  by  initial  pressure.  This  difference  is  plotted  in 
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Figs.  51  and  52  for  two  sets  of  experiments.  On  each  graph,  a  best-fit  curve 
of  the  form: 

Xinv  "  X  "  C1  +  C2  P0  ^  (5'2) 

is  plotted.  This  form  is  anticipated  since  the  boundary-layer  size  is 
affected  by  initial  pressure  or  density  in  that  way  (see  App.  D).  There  is  a 
fair  amount  of  scatter  in  the  data,  but  the  trend  is  quite  clear.  As  the 
initial  pressure  gets  lower,  the  viscous  forces  increase  and  the  difference 
between  the  inviscid  prediction  of  x  and  the  experimentally  measured  value 
increases.  Theoretically,  as  the  initial  pressure  gets  extremely  large,  the 
viscous  forces  should  become  negligible  and  the  solution  should  approach  the 
inviscid  case.  Extrapolating  the  best-fit  curves  to  Pq  =  °°  gives: 


for  e 

-  48° 

X  . 

-  X  =  0.64 

w 

mv 

for  0 

II 

^1 

o 

X. 

-  X  =  -0.95 

w 

mv 

The  fact  that  a  non-zero  value  is  obtained  may  be  due  in  large  part  to 
the  scatter  of  data  involved.  The  scatter  itself  may  be  affected  by  the 
distance  up  the  wedge  that  the  shock  has  progressed  when  the  interferogram  was 
taken  (see  below).  Hornung  and  Taylor  [19]  did  24  experiments,  and  were  able 
to  extrapolate  their  results  more  closely  to  the  inviscid  RR«»>MR  transition 
boundary.  It  is  anticipated  that  if  many  more  experiments  were  done  here,  a 
more  accurate  curve  might  be  drawn,  with  less  scatter  of  data. 

The  influence  of  the  boundary  layer  on  the  triple-point-trajectory  angle 
x  should  De  dependent  on  their  proximity  to  each  other.  As  the  Mach  stem 
grows,  it  gets  further  from  the  boundary  layer.  If  the  boundary  layer  is 
thought  of  as  a  sink,  its  influence  on  x  should  drop  off  with  distance. 
Therefore,  as  the  Mach  stem  grows,  the  boundary-layer  influence  diminishes,  so 


that  the  growth  of  the  Mach  stem,  and  therefore  ;<>  is  inherently  non-linear 
with  time.  The  existence  of  viscosity  eliminates  the  self-similarity  that 
exists  for  the  inviscid  process.  Any  comparison  of  the  triple-point  trajectory 
angle  with  inviscid  predictions  should  be  done  with  respect  to  initial  pressure 
especially  when  the  initial  pressure  and  x  are  quite  small. 


Chapter  6 
CONCLUSIONS 


Many  RR-«-»>MR  transition  boundary  points  in  air  were  determined 
experimentally  by  performing  shock-tube  runs  very  close  to  the  boundary. 
Reflected  wave  angles  were  found  to  be  consistent  with  the  'detachment* 
criterion  for  pseudostationary  flow.  As  anticipated,  initial  low  pressures 
enhance  viscous  effects.  As  the  pressure  is  lowered,  the  transition  boundary 
shifts  downward  on  the  (Ms,  o  )  plane,  giving  rise  to  the  'von  Neumann 
paradox',  in  agreement  with  many  experimenters.  Viscous  effects  have  been 
found  to  be  the  cause  of  the  'von  Neumann  paradox'. 

By  reducing  the  experimental  results  to  a  common  initial  pressure,  a  new 
viscous  RR-*-**MR  transition  boundary  for  air  was  drawn  and  found  to  be  in  good 
agreement  with  the  100  experiments.  The  new  transition  boundary  does  not 
appear  to  be  valid  for  gases  with  low  values  of  y  such  as  Freon-12  and  SF^  at 
high  Mach  number. 

For  a  given  pressure,  the  boundary  layer  has  a  greater  effect  on  the 
transition  boundary  at  lower  Mach  number,  even  though  the  physical  size  of  the 
boundary  layer  is  smaller.  It  is  hypothesized  that  the  degree  of  influence  is 
inversely  proportional  to  flow  deflection  through  the  incident  shock  wave,  due 
to  the  relative  size  of  the  mass  flow  toward  the  wedge  surface  behind  the 
incident  shock  and  the  mass  flow  in  the  boundary  layer. 

The  model  used  to  approximate  the  influence  of  viscous  forces  could  not 
be  substantiated  experimentally,  and  is  assumed  to  be  approximate.  It  did 
however  indicate  that  the  characteristic  length  of  the  flow  is  of  the  order  of 
0.01  mm,  and  therefore  that  the  assumptions  of  laminar  boundary  layer  and 
frozen  flow  as  used  herein  were  appropriate.  It  is  felt  that  correct  and 
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accurate  modelling  of  the  reflection  process  will  probably  come  about  through 
numerical  codes  which  incorporate  the  influence  of  viscosity  through  a 
solution  of  the  Navier-Stokes  equations. 

While  not  directly  investigated,  results  from  MR  experiments  demonstrate 
that  viscosity  also  affects  the  size  of  the  Mach  stem.  Lower  initial  pressure 
reduces  the  size  of  the  Mach  stem,  but  scatter  of  the  limited  data  does  not 
make  for-  a  precise  quantitative  assessment.  This  is  one  area  of  study  which 
should  be  investigated  further. 

In  the  most  general  sense,  viscosity  plays  a  ve» y  significant  role  in 

oblique-shock-wave  reflection.  The  RP-*-*»MR  transition  boundary  is  lowered  a 

significant  amount  (up  to  5),  and  the  Mach  stem  size  in  MR  is  reduced. 

Previous  studies  which  investigated  oblique-shock-wave  reflection  without 

consideration  of  viscous  effects  are  thereby  missing  an  important  factor. 

Definition  of  an  experimental  point  in  terms  of  M  and  ■■  is  insufficient,  for 

j  w 

oblique-shock-wave  reflections  are  a  function  of  M$ ,  u  ,  and  initial  pressure. 
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Figure  2:  Single-Mach  Reflection  (SMR) 


(a)  Schematic  diagram 

(b)  Interferogram  (Ref.  29) 

(■■  =  20 : ,  M  =  3.00,  pQ  =  20.0  kPa,  TQ  =  299  K) 

I  -  incident  shock  wave,  R  -  reflected  shock  wave, 

M  -  Mach  stem,  S  -  slipstream,  T  -  triple  point, 
n  -  wedge  angle 


Figure  3:  Complex-Mach  Reflection  (CMR) 

(a)  Schematic  diagram 

(b)  Interferogram  (Ref.  29) 

(0  =  30' ,  M  -  3.09,  pQ  =  8.0  kPa ,  TQ  =  296  K 

6  -  wedge  angle,  I  -  incident  shock  wave, 

R^  R'  -  first  and  second  reflected  shock  waves, 

M  -  Mach  stem,  S  -  slipstream,  T  -  triple  point,  K 
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Double-Mach  Reflection  (DMR) 

(a)  Schematic  diagram 

(b)  Interferogram  (Ref.  29) 

(  w  =  40  ,  Ms  =  5.50,  pQ  =  1.33  kPa,  TQ  = 

-  wedge  angle,  I  -  incident  shock  wave, 

R,  R'  -  first  and  second  reflected  shock  waves 

M,  M'  -  first  and  second  Mach  stems, 

T,  T'  -  first  and  second  triple  points, 

S,  S'  -  first  and  second  slipstreams 
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Figure  7:  Schematic  diagram  of  regular  reflection  (RR) 


Figure  9:  Schematic  diagram  of  single-Mach  reflection  (SMR) 

I  -  incident  shock  wave,  R  -  reflected  shock  wave, 

M  -  Mach  stem,  S  -  slipstream,  T  -  triple  point, 

9  -  wedge  angle,  X  -  triple-point  trajectory  angle 

w 

<j>  -  incidence  angle,  6  -  deflection  angle 
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Figure  10:  Schematic  diagram  of  complex-Mach  reflection  (CMR) 

I  -  incident  shock  wave,  R  -  reflected  shock  wave, 

M  -  Mach  stem,  S  -  slipstream,  K  -  kink,  T  -  triple  point 
6  -  wedge  angle,  X  -  triple-point  trajectory  angle. 
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Figure  11:  Schematic  diagram  of  double-Mach  reflection  (DMR) 

9  -  wedge  angle,  I  -  incident  shock  wave, 

w 

R,  R'  -  first  and  second  reflected  shock  waves, 

M,  M1  -  first  and  second  Mach  stems, 

T,  T'  -  first  and  second  triple  points, 

S,  S'  -  first  and  second  slipstreams, 

X,  X'  -  first  and  second  triple-point  trajectory  angl 
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Figure  12:  Regular  to  Mach  reflection  transition  boundary  as 
defined  by  two  theoretical  inviscid  criteria 
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Figure  14:  Deviation  of  experimental  results  of  Smith  [8] 

from  the  'detachment'  criterion  RR-*»tMR  boundary. 

Experiments  were  not  performed  in  gaps  between  arrows, 

—  -  .►  lower  limit  of  RR 
- upper  limit  of  MR 
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Figure  15:  Deviation  of  experimental  results  of  Deschambault  [29] 
from  the  'detachment'  criterion  boundary 
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Figure  17:  Effect  of  equilibrium  vibrational  excitation 

of  0^  and  on  the  RR-*»-MR  transition  boundary 
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Boundary-layer  velocity  distribution  in 
moving  and  stationary  reference  frames  [37] 

(a)  Laboratory-f ixed  reference  frame 

(b)  Shock-fixed  reference  frame 

5  -  boundary-layer  height 

x*  -  boundary-layer  displacement  thickness 
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-  based  on  Mirels  [37,  38] 
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Figure  24:  Physical  model  of  the  RR  reflection  process  at  point  P 
in  the  presence  of  a  laminar  boundary  layer  [19] 
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Figure  26:  Predicted  boundary-layer  displacement  thickness  6* 
at  distance  of  1  mm  from  regular  reflection  point  P 
(p0  =  2.0  kPa,  Tq  =  300  K) 


Figure  28:  Effect  of  displacement  wedge  ang 
wave  angle  w'  at  RR-*»MR  transit 
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Figure  30:  Inviscid  solution  for  flow  Mach  number  in  region  1 
at  RR-**-MR  transition  boundary 
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Figure  32:  Experimental  results  for  42^  wedge 
for  high  Pq  and  1 jw  M$ 

Numbers  beside  symbols  are  pn  in  kPa 
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Figure  34:  Experimental  results  for  47  wedge 

for  low  pn  and  M 
U  s 

Numbers  beside  symbols  are  Pq  in  kPa 
O  RR 
□  SMR 
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Figure  35b:  Experimental  results  for  48“  wedge 
for  low  Pq  and  M$ 

Numbers  beside  symbols  are  Pg  in  kPa 
O  RR 
□  SMR 
A  CMR 
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Fiqure  36: 


Experimental  results  for  48  wedge 
for  low  pQ  and  high  M$ 

Numbers  beside  symbols  are  Pq  in  kPa 
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Figure  38:  Experimental  results  for  45  wedge 
for  low  Pq  and  high 

Numbers  beside  symbols  are  Pq  in  kPa 
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Figure  41b: 


Shadowgraph  of  complex-Mach  reflection  (CMR) 
Uoseup  of  triple  point  T 
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Analytic  prediction  of 
triple-point  trajectory 
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Figure  43: 


Comparison  of  experimentally  determined  transition 
boundary  points  with  inviscid  RR«**-MR  boundary 
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Figure  46:  Displacement  wedge  angle  0^  at  RR-**»MR  transition 

boundary  for  Pq  =  2.0  kPa  as  determined  by  reduction 
of  experimental  results  at  various  initial  pressures 
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Figure  48:  Characteristic  length  of  boundary  layer  as  determined 
by  comparison  of  predicted  boundary-layer  size  and 
experimental  results 
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Ratio  of  boundary-layer  size  parameter  si  to 
experimentally  determined  displacement  wedge  angle 
vs  incident  shock  flow  deflection  angle  0^ 


Figure  50a:  Predicted  viscous  RR-*^MR  transition  boundary 
in  air  for  initial  pressure  pn  =  5.0  kPa 
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Figure  50b:  Predicted  viscous  RP-*-*41R  transition  boundary 
in  air  for  initial  pressure  pn  =  5.0  kPa 
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Figure  52:  Difference  between  analytic  triple-point-trajectory 
angle  X_nv  and  experimentally  measured  trajectory 
angle  X  for  0  =  48 


Appendix  A 


SOL UTIQN  OF  FLOW  FOR  RR  AND  M R 

The  flow  at  an  intersection  of  shock  waves  is  readily  calculated  by 
applying  the  equations  of  motion  across  each  of  the  waves,  and  using  the 
proper  boundary  conditions.  The  frame  of  reference  is  fixed  to  the  confluence 
point;  triple  point  T  in  MR  and  reflection  point  P  in  RR.  The  equations  of 
motion  across  a  shock  wave  are  (Figs.  A-l  and  A-2): 

Continuity  of  tangential  velocity: 

U-  cos(4>.)  =  U.  cos(4>.-0,)  (A.l) 

1  1  J  J 

Conservation  of  mass: 

p.  Ui  sin(4>i)  =  p.  Uj  sin(<t>i-0.)  (A. 2) 

Momentum  equation: 

p.  +  Pi  Ui2  si n2(4>i )  =  Pj  +  Pj  CL2  sin2(4>i~Gj )  (A. 3) 

Conservation  of  energy: 

hi  +  |  s  i  n 2  ( <p  i )  =  hj  +  l-  V2  sin2^-^ )  (A. 4) 

where  subscripts  i  and  j  refer  to  the  regions  upstream  and  downstream  of  the 
shock  wave  respectively. 

The  two  necessary  independent  thermodynamic  variables  usually  specified 
are  p  and  T,  thus  giving  p  and  h  since  the  equations  of  state  are  given  by: 

,■  =  p(p,T)  (A. 5) 

h  =  h(p,T)  (A. 6) 

If  the  gas  is  frozen,  a  closed-form  solution  exists. 


A1 


For  both  RR  and  MR  the  following  substitutions  apply: 

Incident  shock  wave  I:  i  =  0  ;  j  =  1  (A. 7) 

Reflected  shock  wave  R:  i  =  1  ;  j  =  2  (A. 8) 

Mach  stem  M:  i  =  0  ,  j  =  3  (A. 9) 

(except  4>i  =  <t3) 

For  clarity,  these  are  illustrated  in  Fig.  A-l  for  the  case  of  MR,  and 
Fig.  A-2  for  RR. 

The  variables  most  often  used  to  define  the  case  under  consideration  are 

M  and  0  .  The  incident  flow  velocity  and  direction  are  calculated  by 
s  w 

attaching  the  frame  of  reference  to  point  T  in  MR,  or  P  in  RR,  and  applying 
the  gasdynamic  equations  and  geometric  relationships. 

Solution  of  Flow  at  T  in  MR 

The  boundary  conditions  are: 

P2  =  P3  (A. 10) 

9j  -  e2  =  e3  (a. ii) 

An  additional  constraint  must  be  specified  to  solve  the  equations.  It  is 
assumed  that  at  T,  the  Mach  stem  is  perpendicular  to  the  wedge  surface. 

One  method  of  solution  is: 

1.  Guess  X  and  w' ,  and  calculate  the  flow  in  region  0, 
relative  to  T. 

2.  Calculate  the  flow  in  region  1  by  applying  the  equations  of 
motion  across  I. 

3.  Calculate  the  flow  in  region  2  by  applying  the  equations  of 
motion  across  R. 


A2 


4.  Calculate  the  flow  in  region  3  by  applying  the  equations  of 
motion  across  M. 

5.  Compare  the  flow  direction  and  pressure  in  regions  2  and  3, 
and  adjust  X  and  uj  1  to  try  and  achieve  boundary  conditions 
A. 10  and  A. 11. 

6.  Repeat  steps  2  through  5  until  convergence  is  achieved. 

It  should  be  noted  that  in  some  areas  of  the  (Ms,  ©w)  plane,  two 
solutions  for  MR  can  exist  [44],  However,  the  one  with  lower  value  of  X  and 
uj  1  is  the  only  one  which  occurs  in  practice. 

Sol ution  of  Flow  at  P  in  RR 

This  configuration  is  illustrated  in  Fig.  A-2.  The  only  boundary 
condition  is: 

e1  -  02  =  0  (A. 12) 

The  method  of  solution  is  very  similar  to  that  for  T  in  MR. 

1.  Guess  u',  and  calculate  the  flow  in  region  0,  relative  to  P. 

2.  Calculate  the  flow  in  region  1  by  applying  the  equations  of 
motion  across  I. 

3.  Calculate  the  flow  in  region  2  by  applying  the  equations  of 
motion  across  R. 

4.  Find  the  flow  direction  in  region  2  and  adjust  w1  to  try  and 
achieve  boundary  condition  A. 12. 

5.  Repeat  steps  3  to  4  until  convergence  is  achieved. 

Solutions  for  RR  and  MR  over  a  wide  range  of  Mach  numbers  and  wedge 
angles  were  tabulated  and  plotted  by  Hu  and  Shirouzu  [27]. 


APPENDIX  B 


I 


Computer  Program  to  find  the 


RR-^-MR  Transition  Boundary  in  Frozen  Air 


*»*»*'»«>**' *4  -»**»««'»*  *  **•«  *  **■»■■»  4-iH*  «**«•»  #«"»  ■»  4  ***  * 

*  THIS  PROGRAM  SOLVES  FOR  THE  RR-MR  TRANSITION  BOUNDARY  FOR  FROZEN  H 

4  AIR  IT  CAN  BE  EASILY  MODIFIED  FOR  OTHER  FROZEN  GASES.  * 

*  OPTIONS  ARE  FOR  THE  DETACHMENT,  SONIC,  AND  MECHANICAL  EQUILIBRIUM  * 

*  CRITERIA  -i 

*  THE  SHORT  AND  SIMPLE  DESIGN  ALLOWS  IT  TO  BE  TRANSLATED  FOR  A  * 

*  MICROCOMPUTER,  AND  EXECUTION  TIME  IS  MINIMAL  * 

4  **■»•■*  •»■»•>)  rt  ****  4-tm-U  *■»■■»*  *«#«  4  ■»  ##■###■»■#■»#  4 -*4  -tH 

REAL  MS,  MO,  Ml  ,  M2 
COMMON  G,  R 
PI =3  141592654 
RADDEG-180  O/PI 


THERMODYNAMIC  GAS  PROPERTIES 

R  =8320  1/2S  967 
G- 1  400 


OPTIONS 


I  OPT 
IOPT 
IQPT 


DETACHMENT  CRITERION 
SONIC  CRITERION 

MECHANICAL  EQUILIBRIUM  CRITERION 


READ*  5.  ■»  >  NCA5E 
WRITE  (6,  100) 

DO  10  1=1, NC ASE 

THETAD  IS  THE  DISPLACEMENT  WEDGE  ANGLE  CAUSED  BY  THE 
BOUNDARY-LAYER  BEHIND  THE  INCIDENT  SHOCK  WAVE 
IT  APPLIES  ONLY  TO  THE  SONIC  AND  DETACHMENT  CRITERIA 
TRANSITION  BOUNDARIES 

READ*  5-*  )  MS,  THETAD,  IOPT 
IF* IOPT  EQ  3)  THE  T  AD  =  0  0 
T HET A D= THETAD /RAD DEG 

I F ( IOPT  LE  2)  CALL  DE TSON ( MS,  THET AW,  THETAD, PH  1 1 ,  I  OPT > 
IF*  IOPT  EQ  3)  CALL  DETSON  (  MS,  THET  AW.  0  0,  PHI  1,1) 

IF*  IOPT  EQ  3)  CALL  MECHEQ  *  MS,  THET  AW,  F  H 1 1 ,  I  OPT  ) 

CALCULATE  ALL  THE  THERMODYNAM IC  STATES 

PH10=F'l/2  -THETAW 
MO=MS  /  S  I  N  (  PH  I  0 

CALL  OSWAVEIPMIO.  MO,  1  ,1  ,  Ml ,  P  10 .  T  10 .  THET  A 1  ) 


noon  non  n  n  n  n  n  n  n  n  n  non 


CALL  OSWAVECPHIl,  Ml,  P  10,  T  10,  M2,  P20,  T20,  THETA2  ) 

0MEGAP=PH 1 1-THETA  1 

SWITCH  FROM  RADIANS  TO  DEGREES  AND  OUTPUT  THE  DATA 

THETA W=THETAW*RADDEG 

PHIO=PHIO*RADDEG 

PHI1=PHI1*RADDEG 

THETA  1 =THETA1 *PADDEG 

THETA2=THETA2*RADDEG 

THETAD=THETAD*RADDEG 

OMEGAP  =0MEGAP  *RAD DEG 

WR ITE  (  6,  200)  MS,  THETAW,  THETAD,  lOPT.MO,  PHIO,  THETA  1 ,  Ml,  P10,  T10, 
1PHI1,  THETA2,  M2,  P20,  T20,  OMEGAP 
10  CONTINUE 
STOP 


100  FORMAT  (1H1,  ///,  T7,  'MS  ',  T13,  'THETAW',  T21,  'THETAD',  T30,  'IOPT  ', 

1T41,  'MO',  T 48,  'PHIO',  T55,  'THETA1  ',  T65,  'Ml  ',  T72,  'P1/P0',  TSO,  'Tl/TO', 
2T88,  '  PHI  1  ' ,  T95,  'THETA2 T  105,  'M2',TU2,  'P2/P0  ' ,  T1  20,  'T2/T0', 

3T127,  'OMEGAP ' ,  ///  ) 

200  FORMAT  ( 1H  ,  T5,  Fh.  3,  T13,  F6.  2,  T21,  F6.  2,  T32,  II,  T39,  F6.  3,  T47,  F6  2, 

1 T55,  F 6  2,  T 63,  F6  3,  T71 ,  F6  2,  T79,  F6.  2,  T97,  F6  2,  T95.  F6  2,  T103,  F6  2, 

2T 111,  F6  2,  T119.F6  2,  T127,  F6  2> 

END 

*#***»*»»#*##*#+****#»***»#***m#m#m*t*«M*»*****m»*#*#*##«* 
FUNCTION  PHI M ( Ml ) 

THIS  FUNCTION  CALCULATES  INCIDENCE  ANGLE  FOR  MAXIMUM  DEFLECTION 
THROUGH  THE  SHOCK  WAVE 

REAL  Ml 
COMMON  G,  R 

C  1  =  (  G  + 1  )  /4.  *  (Ml**2)-1. 

C2=SQRT((G+1  )*(1.+(G-1  ) /2  *  C  Ml  *  *2  )  +  (  G+  1  > / 1 6  * < M 1**4 > ) > 

PHIM=ASIN ( SORT  < 1  /G/ ( M 1**2 ) * <  C 1 +C 2 > )  > 

RETURN 

END 

*##**#*#tHHHH*»1HHMHHt**'******#»*«*iHHt»**»#»*****«  -tMHMt  •»*•«•■»■»*******  ** 


FUNCT ION  PHIS (Ml ) 

THIS  FUNCTION  CALCULATES  THE  INCIDENCE  ANGLE  FOR  SONIC  FLOW 
BEHIND  THE  SHOCK 

REAL  Ml 
COMMON  G .  R 

C  1  =( G+ 1  )  / 4  *  <M1**2>- <3.  -G>/4 

C  2=SQR  T  f  f  G  +  l  >*((9  +G  >  / 16  -<3  -G)/8  *  <  Ml  *  *2  >  +  <  G+ 1  )/16  *(M1**4)>) 

PHIS -AS IN ( SQR  T  < 1  /G/ ( Ml**2 )* ( C 1+C 2 ))  ) 

RETURN 

END 


B  2 


non  n  on  nooonnoooo  non 


C 

r 

C 


fr**#*#^***-#-****####*#*#*#  ■»*•«■■»  *■#■#■»  #«##  «“»*"*  *■*#■■»  ■»# 


subroutine  oswaveiPhi  i,  mi,  pi,  ti,  m2,  P2,  T2,  theta > 

FLOW  PROPERTIES  ARE  CALCULATED  THROUGH  A  SHOCK  WAVE 

REAL  Ml,  M2,  MIN,  M2N 
COMMON  G,  R 
M1N=M1*3IN<PHI 1 ) 

M2N=SQRT(  <2  +  'G-1  >•*<  M1N**2)  )  /  <2.  *G*<  M1N**2)-G  +  1  )) 

P2=P1*(2.  *G/(G-*-l.  >*(M1N**2)-(G~1.  >/(G  +  l  )) 

T2=T1*< 'l.+IG-l  ) /2  *  <M1N**2)  ) /( 1 .  +(G-1.  )/2  * <M2N**2)  )  ) 
Cl  =  i5+1  ) /2 .  * CM1**2>/ ( (M1N**2 )-l  )-l 

THETA-ATAN(1.  /C 1 / TAN ( PHI  1  )  > 

M2=M2N/SIN(PHI 1-THETA ) 

RETURN 

END 


*  ****■»■■»•&■■«■•»*■!»■ -a  lHHHHHHHt  *■  a***  a  *-)nt  «»*«*«*###»***«« 


SUBROUTINE  DETSON(MS,  THETAW,  THETAD,  PH  1 1,  I  OPT) 

THIS  SUBROUTINE  SOLVES  FOR  EITHER  THE  DETACHMENT  OR  SONIC  CRITERION 
THE  FOLLOWING  LOGIC  IS  USED: 

1  SET  THE  MAXIMUM  AND  MINIMUM  POSSIBLE  VALUES  FOR  PHIO 

2  ITERATE  ON  PHIO  BY  BISECTING  THE  RANGE  OF  PHIO  UNTIL 

THE  LIMIT  OF  RR  IS  REACHED,  WHETHER  IT  BE  THE  LIMIT  IMPOSED 
BY  THE  DETACHMENT  OR  SONIC  CRITERIA 


REAL  MS,  MO,  Ml  ,  M2 
COMMON  G.  R 
EXTERNAL  PHIM.PHIS 
P I =3  141592654 
N--Q 

MINIMUM  POSSIBLE  VALUE  FOR  PHIO 
P  H IQA =0  0 

MAXIMUM  POSSIBLE  VALUE  FOR  PHIO 
PH IOB  =P I / 2 

LOOK  FOR  BOUNDARY 

10  PHIGT  =  (PHI  OA  +  PHI0B)/2 
MO=MS /SIN ( PHI OT ) 

CALL  QSWAVE ( P  H IOT  ,  MO ,  1  ,  1  ,  M 1  ,  P 1 ,  T 1 ,  THETA  1 ) 
I F ( M 1  GE  1  )  GO  TO  20 
C  RR  NOT  POSSIBLE  -  PHIO  TOO  HIGH 
F  UNC  T  =  1  0 
GO  TO  30 

20  PHI 1=PHIM (Ml  ) 

IFOOPT  EQ  2)  PHI  1=PHI3(M1  ) 

CALL  QSWAVE(PHI1,  Ml,  1  ,1  ,  M2,  P2,  T2,  THETA2  ) 
FUNCT  =THETA1 - THET  A2-THETAD 


BO 


on  nr  j 


N=N+ 1 

r F (N  GT  20)  GO  TO  50 
C 

MAKE  A  NEW  GUESS  FOR  PHIO  USING  THE  BISECTION  METHOD 

30  IF(FUNCT  GT  0  0)  THEN 
P  H IOB  =PH I OT 
ELSE 

PHIOA^PHIOT 
ENDIF 
GO  TO  10 

50  THETAW=P 1/2  -PHIOB 
RETURN 
END 

##*# •'Mi  »1HHt*#4****^HHHt#*#**#»#m#*#*#*#*#»**#*#*#***** 

C 

SUBROUTINE  MECHEQ  (  MS,  THETAW,  PHI  1 ,  IOPT) 

C 

C  THIS  SUBROUTINE  SOLVES  THE  MECHANICAL  EQUILIBRIUM  CRITERION 
C  THE  FOLLOWING  LOGIC  IS  USED: 

C 

C  1  SET  THE  MINIMUM  AND  MAXIMUM  POSSIBLE  VALUES  FOR  PHIO. 

C  2  CHECK  TO  SEE  IF  A  SOLUTION  IS  POSSIBLE  BY  LOOKING 

C  AT  THE  DETACHMENT  CRITERION  SOLUTION. 

C  THIS  POINT  CORRESPONDS  TO  THE  MAXIMUM  POSSIBLE  VALUE 

C  FOR  PHIO  AT  THE  MECHANICAL  EQUILIBRIUM  BOUNDARY. 

C  3  USE  THE  BISECTION  METHOD  TO  GET  A  FIRST  GUESS  FOR 

C  PHIO  AT  THE  MECHANICAL  EQUILIBRIUM  BOUNDARY. 

C  4  CALCULATE  THE  FLOW  IN  REGIONS  1  AND  3. 

C  AND  FLOW  DEFLECTION  THROUGH  THE  INCIDENT  SHOCK 

C  5  ITERATE  ON  PHI1  UNTIL  THE  DEFLECTION  THROUGH  THE  REFLECTED 

C  SHOCK  IS  EQUAL  AND  OPPOSITE  TO  THAT  THROUGH  THE  INCIDENT 

C  SHOCK  USE  A  BISECTION  METHOD 

C  6  COMPARE  THE  PRESSURES  IN  REGION  2  AND  3. 

C  7  USING  THE  BISECTION  METHOD,  GET  A  NEW  VALUE  FOR  PHIO 

C  AND  REPEAT  STEPS  4  THROUGH  6  UNTIL  CONVERGENCE  IS 

C  ACHIEVED 

C 

REAL  MS-  MO,  Ml  ■  M2,  M3 
COMMON  G,  R 
EXTERNAL  PHIM 
P I =3  141592654 

N-0 

C  MINIMUM  POSSIBLE  VALUE  FOR  PHIO 

PH  IOA  -~0  0 

C  MAXIMUM  POSSIBLE  VALUE  FOR  PHIO,  THE  DETACHMENT  CRITERION  BOUNDARY 
P H IOB  =P 1/2  -THETAW 
C 

C  CHECK  TO  SEE.  IF  A  MECHANICAL  EQUILIBRIUM  SOLUTION  EXISTS 
C 

MO=MS / S I N ( PH  I  OB  > 

CALL  OSWAVE' PHIOB, MO,  1  ,1  , M 1 ,  P 1 ,  T 1 ,  T HETA 1 > 

CALL  OSWAVE I  / 2  ,  MO,  1  ,1  ,  M3 ,  P3,  T3,  DUMMY  ) 

PHI 1A=PHIM(M1 ) 

CALL  GSW4VE(PHI  1  A,  Ml  ,  PI,  T  1 ,  M2,  P2,  T2,  THETA2) 


non  n  n  n  n 


ir-PJ  l.£  PZ  >  GO  TQ  10 

WRITE  6.  *  )  NO  MECHANICAL  EQUILIBRIUM  SOLUTION  EXISTS  AT', 

1  MS  =  ■  MS 

i Of  r  =  i 
RETURN 
C 

C  LOO*  FOR  THE  BOUNDARY 

c 

10  PHITT-i PH  I OA+PHJOB 1 /2 
M O =MS / S I N ' PH  I U T ) 

CALL.  OSWAUEvPmIOT  ,  MO,  1  ,1  ,  Ml  ,  P 1 ,  7 1 ,  THETA  i  ) 

CAl.L.  OSvYAUE'RI  2  .  MO,  1  ,1  ,  M3,  P3,  T3>  THETA3) 

C 

C  SOLVE  Pnc  -‘HU  T  J  ACHIEVE  NO  NET  DEFLECTION  -  USE  BISECTION  METHOD 
C  TO  ITERATE 

C 

N2  =  u 

C  MAXIMUM  POSSIBLE  VALUE  OF  PHI1 
PHI i A=PHIM( Ml > 

C  MINIMUM  POSSIBLE  VALUE  OF  PH  II 
PH  I  1 B=ASIN (  I  /Ml >  +  .000001 
20  PHI  1T=<  PHI1A+PHI1B)  /2. 

CALL  OSWAVE  (  PH  I  1  T,  M  1 ,  P 1 ,  T 1 ,  M2,  P2,  T2,  THETA2T  > 

IF  <  ABS  >  THETA2T -THETA!)  LT.  1.7E-05)  GO  TO  30 
IF ( N2  GT  20 )  GO  TO  40 
N2=N2+1 

IF ( THET  A2T  GT  THETA  1 >  THEN 
PHI 1A=PHI IT 
ELSE 

PHI 1B=PHI IT 
END  IF 
GO  TO  20 

CALCULATE  A  NEW  VALUE  FOR  PH  10 
USE  THE  BISECTION  METHOD  TO  ITERATE 

30  CONTINUE 

FUNCT  =P3~P2 
N  =N+ 1 

IF ( N .  GT  20)  GO  TO  60 
I F (FUNCT  GT  0.  0)  THEN 
PHIOA=PHI OT 
ELSE 

PHIOB=PHIOT 
ENDIF 
GO  TO  10 

ERROR  MESSAGES 

40  WRITE (6,  *)  'NO  CONVERGENCE  -  INNER  LOOP  OF  MECHEQ  ' 

STOP 

60  THETAW=PI/2  -PHIOB 
PHI 1=PHI 1 T 
RETURN 
END 


B  5 


vv 


It  is 


The  values  of  and  from  two  sources  are  listed  in  Table  C-l. 
important  to  note  that  the  relaxation  tire  is  inversely  proportional  to 
pressure,  and  much  less  dependent  on  temperature. 

The  relaxation  lengths  behind  a  normal  shock  wave  in  0^  and  are 

plotted  in  Figs.  C-l  and  C-2.  A  pressure  of  2.0  kPa  was  chosen  as  a  reference 
value  for  the  plots,  but  correction  for  any  other  pressure  can  be  made  using 
Eq.  C.5. 

One  might  expect  that  the  relaxation  lengths  for  oxygen  and  nitrogen  will 
differ  when  the  two  are  mixed  as  in  air.  Since  oxygen  is  affected  before 
nitrogen,  its  vibrational  mode  is  not  affected  much  by  collision  with 
vibrating  nitrogen  molecules.  As  such,  the  relaxation  length  of  oxygen  is 
affected  very  little  by  the  presence  of  nitrogen  [45].  This  is  not  true  for 
nitrogen  however.  When  nitrogen  begins  to  become  excited,  the  oxygen 
molecules  are  already  vibrating,  and  their  collision  with  the  nitrogen 
molecules  does  change  the  relaxation  length  of  nitrogen.  Even  with  this 
change,  the  relaxation  length  of  oxygen  is  still  much  shorter  than  that  of 
ni trogen. 


Table  C-l 


Empirical  Constants  in  the 
Vibrational-Relaxation-Time  Equation  (C.5) 


Gas 

C1 

(atm-ysec) 

(Deg!  K1/3 ) 

Temp.  Range 
(Deg.  K) 

Ref. 

°2 

2.92  x  10’4 

126 

1000  -  3700 

46 

°2 

5.42  x  10"5 

143.4 

800  -  3200 

47 

n2 

6.22  x  10'5 

202 

3000  -  4500 

46 

n2 

7.12  x  10“3 

124.1 

800  -  6000 

47 

C4 


Hill 


Ref.  46 


Ref.  47 


INCIDENT-SHOCK-WRVE  MRCH  NUMBER  M 


Figure  C-l:  Vibrational  relaxation  length  of  0^  behind  a  norma1 
shock  wave  with  Tq  =  300  K  and  =  2.0  kPa 
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Figure  C-2:  Vibrational  relaxation  length  of  N2  behind  a  normal 
shock  wave  with  Tn  =  300  K  and  pn  =  2.0  kPa 


Appendix  D 


LAMINAR-BOUNDARY-LAYER  EQUATIONS 

This  appendix  outlines  the  basic  technique  used  in  calculating  the 
boundary-layer  profile  behind  a  moving  shock  wave.  This  method  of  solution  is 
based  on  work  done  by  Mirels  [37,38].  Recent,  more  precise  methods  of  solution 
exist  [48,  49]  but  were  not  used  due  to  complexity.  The  reference  frame  is 
attached  to  the  shock  wave,  as  shown  in  Fig.  D-l,  so  that  the  flow  (U  -  U^)  is 
steady,  and  the  wall  moves  with  a  velocity  U  . 

The  equations  which  govern  the  flow  in  the  boundary  layer  are  the  same  as 
those  for  the  external  flow,  continuity,  momentum,  energy  and  an  equation  of 
state.  The  boundary-layer  equations  include  terms  for  viscosity  u  and  heat 
conduction  k,  and  in  this  solution  it  is  assumed  that  no  pressure  gradient 
exists.  For  a  more  general  explanation  of  the  equations  and  their  physical 
meaning,  see  Schlichting  [50]. 


The  boundary  layer  equations  are: 


?  ( r.  u )  +  I*  ( c  v )  =  0 

ax  ay 

Continuity 

(D.l) 

U  5U  +  V  jli  .  1  3  / m  _3U  \ 

ax  ay  "  :  ay  \  ay  / 

Momentum 

(D.2) 

Energy 

(D.3) 

p  =  CRT 

Equation  of  State 

(D.4) 

with  boundary  conditions: 

u(x,0)  =  uw 

(D.5a) 

u(x,«)  =  ue 

(D.5b) 

< 

X 

O 

II 

o 

(D.5c ) 

T(x,0)  =  Tw 

(D.5d) 

T(x.°")  =  T& 

(D.5e) 

D1 


Since  the  flow  is  continuous,  a  stream  function  ij;  exists  such  that 


i 


£V 

?y 


eX 


£U 

pr 

£V 

p~ 


Introduce  a  similarity  parameter 


=  pT~  ft 
"V  2xvr  J  T 


dy 


and  rewrite  the  stream  function 


:^/2uexvr  f ( n) 


It  should  be  noted  that 


=  f* 


T_ 

T. 


(f  +  2><f' t ) 


(D.6a) 

(D.6b) 


(D.7) 


(D.8) 


(D.9a) 

(D.9b) 


As  a  simplification,  it  is  assumed  that  the  viscosity  and  the  thermal 
conductivity  are  directly  proportional  to  the  temperature. 


p  - 


=  T 


(D.lOa) 


k  =  T 


(D. 10b) 


The  choice  of  reference  temperature  for  viscosity  and  thermal  conductivity  is 
somewhat  arbitrary  at  this  point.  A  wise  choice  of  reference  temperature  may 
minimize  error  caused  by  this  simplification  [38],  and  will  be  discussed  later. 


When  equations  D.8,  D.9,  and  D.10  are  substituted  into  the  momentum 
equation,  it  becomes 

f"  +  ff'"  =0  (D.  11) 

with  boundary  conditions 

f  (0)  =  0  (D.lla) 

f 1 (0)  =  uw/ue  (D. lib) 

f '(«)  =  1  (D. 11c) 

Substitution  of  equations  D.9,  D.9,  and  D.10  into  the  energy  equation 
yields: 


T' 1  +  ofT'  =  -o  ( A— 1 )  Me2  (f ")2 
assuming  the  Prandtl  number  a  is  constant. 

From  equation  D-12,  the  temperature  distribution  can  be  expressed  as 


(D.12) 


r =  1  + 

1  e 


Mg  r(n)  + 


where  r(n)  satisfies 


r"  +  ofr'  =  - 


_ 

(H 


(f")2 


(D.13) 


(D. 14) 


s(n)  satisfies 


r  *  ( 0 )  =  0 

(D. 14a ) 

r  (°°)  =  0 

(D. 14b) 

+  of s’  =0 

(D.15) 

s  ( 0 )  =  1 

(D.15a) 

s(°°)  =  0 

(D. 15b) 

>  *J-  - 

The  temperature  distribution  is  in  fact  a  combination  of  2  solutions;  the 

solution  for  the  case  of  an  insulated  wall,  plus  one  for  addition  of  heat 

transfer  through  the  wall.  T  •  is  the  temperature  which  would  occur  at  the 

w,  1 

wall  if  it  were  insulated,  and  r(0)  is  a  recovery  factor  for  the  recovery 
temperature  T^  at  the  wall,  which  is  less  than  the  adiabatic  total  temperature. 

To  get  y  in  terms  of  n,  equations  D.7  and  D.13  are  combined,  giving 


y  = 


I2xv  T 
r  e 


ue  ^  r 


,  +  111 
2 


6 -)']■/:* 

J  o  o 


S  dr,  (D.  17) 


Another  parameter  of  interest  is 
T 


**  =  ?x  — 
T 

r 


r  -i 

1  im  (r-f )  +  Izl 

(f  -  >)  He 

2xu* 

e 

r.  co 

\  e  / 

> 

7 


r  dn 


rw  ~  Tw,i  f 


n 

s  dn 


(D.18) 


The  distribution  of  6  with  axial  distance  is  of  the  form: 


6*  =  Cj 


(D. 19) 


The  angle  at  which  flow  enters  the  boundary  layer  is  determined  by  the  rate 
at  which  *  grows.  The  boundary-layer  flow  entry  angle  is  defined  by: 


Solution  of  equations  D.ll  to  D.llc  and  D.14  to  D.15b  gives  a 
non-dimensional  profile  of  the  boundary  layer.  The  rest  of  the  equations  can 
then  be  applied  to  get  velocity  and  temperature  profiles. 

The  wall  temperature  used  is  that  of  the  wall  before  the  shock  passes. 
This  was  shown  to  be  quite  valid  [38]  as  the  conductivity  of  the  wall  is  much 
higher  than  that  of  the  air,  and  the  wall  acts  as  a  large  heat  sink. 

Viscosity  and  heat  conduction  do  not  vary  linearly  with  temperature  and 
the  assumption  that  they  do  will  only  be  approximate.  To  minimize  the  effects 
of  this  approximation,  a  specific  reference  temperature  is  used.  The  reference 
temperature 

Tr  =  0.5  (Tw  +  Te)  +  0.22  (TwJ  -  Tg)  (D.19) 
will  minimize  the  error  in  heat  transfer  and  shear  stress  at  the  wall  [38], 

The  method  of  solution  for  the  non-dimensional  profiles  is  as  follows: 

1.  Solve  for  the  distribution  of  f 

(a)  Set  boundary  conditions  D.lla  and  D.llb  for  f(0)  and 
f'(0) 

(b)  Guess  f''(0)  and  use  a  Runge-Kutta  method  to  integrate 

equation  D.ll  from  n=0  to  n=6  (n=6  is  a  close 

approximation  to  n  =  00 ) 

(c)  Compare  f(6)  to  boundary  condition  D.llc 


(d)  Repeat  steps  b  and  c  until  convergence  is  achieved 


2.  Solve  for  the  distribution  of  r 


(a)  Set  f(0),  f'(0),  f ' 1 ( 0 )  obtained  in  step  1 


(b)  Set  boundary  condition  D.14a  for  r'(0) 


(c)  Guess  r(0)  and  use  a  Runge-Kutta  method  to  integrate 
equations  D.ll  and  D.14  from  n  =  0  to  n  =  6. 


(d)  Compare  r(6)  to  boundary  condition  D.14b 


(e)  Repeat  steps  c  and  d  until  convergence  is  achieved. 


3.  Solve  for  the  distribution  of  s 


(a)  Set  f( 0),  f  1 ( 0 ) ,  f ' ' ( 0 )  obtained  in  step  1 


(b)  Set  boundary  condition  D.15a  for  s(0) 


(c)  Guess  s ' (0)  and  use  a  Runge-Kutta  method  to  integrate 
equations  D.ll  and  D.15  from  n  =  0  to  n  =  6. 


(d)  Compare  s ( 6 )  to  boundary  condition  D . 15b 


(e)  Repeat  steps  c  and  d  until  convergence  is  achieved. 


The  complete  Fortran  program  for  the  solution  of  the  boundary  layer 
equations  is  listed  in  Appendix  E. 


//v.  f.  r.  s.  -•  -■ 
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APPENDIX  E 


Computer  Program  for  the  Solution  of  the 
Laminar— Boundary-Lager  Equati ons 


SUBROUTINE  BLAYER  (PE,  TE,  UE,  TW,  UW,  REYN,  DELTA,  DELSTAR,  THETAD, 

1  I  OUT  • 

C 

C  THIS  SUBROU’ INE  SOLVES  FOR  THE  BOUNDARY  LAYER  BEHIND  A  PLANE 
NORMAL  SHOCK  WAVE  MOVING  INTO  STILL  AIR 

IT  IS  BASED  ON  NACA  TN  3401  AND  NACA  TN  3712  BY  H  MIRELS 

REAL.  ME,  NOW,  NUREF  ,  A<  301,  6  ) 

R  EAL  Y  o  >  ,  UIM  1  00  >  BLR  (  30  1 , 7 ) 

INTEGER  I WK ( 5  > 

COMMON  ■  C  /  PR  Cl. 

EXTERNAL  FCN,  FCNj 

INI  T  IAL,I  ZE  DATA 

PP= 0  72 
C  1  =0  0 

R  A IR  =  28  7  074 
G  A  { R  =  1  4 

PW=PE 

CLENGTH  IS  AN  ARBITRARILY  CHOSEN  CHARACTERISTIC  LENGTH 
C  LP NGTH= 1  OE— 0  3 

CALCULATE  SOME  PLOW  PROPERTIES 

AE-SGRT  GA  IR-a-RA  I  R  *TE  > 

MC=UE /AE 
R HOE  =  PE  'PAIR  /  TE 
R  HOW  =PW/RAIR/TW 
N  U W~  V  I  SC  U  (  TW  )  /  RHO  W 
UWUE  =  UW/:JE 

*«««*«*■*  SOLUTION  OF  BOUNDARY  LAYER  DIFFERENTIAL  EQUATIONS 

Yfl)  =  F 

Y  <  2  >  -  r 

Y  ( 3  )  =  F'1 

Y  ( 4  •  --  R  OF  S  DEPENDING  ON  WHICH  IS  BEING  SOLVED  FOR 
Y (5,  -  P  OP  S'  DEPENDING  ON  WHICH  IS  BEING  SOLVED  FOR 

*■**#■•»  «  w  »  «  »  4 

*  SOLVE  FOR  THE  I)  I STR  I  BUT  1  ON  OF  F  * 


nonon  nnn  n  non  oonoo  no  non 


INITIALIZE  DATA  FOR  IMSL  SUBROUTINE  'DGEAR ' 

'DGEAR '  INTEGRATES  THE  BOUNDARY  LAYER  EQUATIONS  USING  RUNGE-KUTTA 
H.  N,  TOl,  METH,  AND  MITER  ARE  PARAMETERS  FOR  'DGEAR  ' 


N  =  5 

T  QL=  OOOOOl 
METH=2 
MI  TER  =3 

SET  THE  INITIAL  CONDITIONS 
FPPZ  IS  F " ( 0  > 

FPFZ=-0  1  24425*  (  UWUE**2)  -0.  77505*UWUE  +  1.  0287 
10  X  =0  0 

y  ( i  >=o.  o 

Y ( 2 ) =  UWUE 
Y ( 3 ) =FPP  Z 
Y  (  4  )  =  0 .  0 
Y <  5>=0  0 
H=  OOOOOl 
I NDE  X  =  1 

SOLVE  FOR  F-PRIME  AT  INFINITY  (FPI) 

DO  20  1=1 , 300 
XEND=. 02* I 

CALL  DGEAR  (N/  FCN,  FCNJ,  X,  H,  Y,  XEND,  TOL,  METH,  MITER,  INDEX,  IWK,  WK,  IER  ) 
20  IFvIER  NE.O)  WRITE<6. *>  'IER  IS  ' , IER 

FPPZNEW=FPPZ 
FP INEW=Y ( 2 ) 

I  F  (  IOUT  EQ.  2  )  WRITE<6,*>  'FPINEW  IS  '  ,  FP  I  NEW,  FPPZ 
L=L+1 

I F  ( L  GT  15)  STOP 
I F  <  L  GT  1  )  GO  TO  30 
FPP  Z  =  FPP  ZNEW+O.  05 
GO  TO  40 

CALCULATE  A  NEW  F" (0)  —  FPPZ 

30  I  F  ( AB  S  ( FP  I NEW-1  0  ) .  LT.  1.  0E-05  )  GO  TO  100 

FPPZ  =  FPPZNEW+ ( 1.  O-FPINEW) / (FP I NEW-FP I  OLD > * (FP PZNEW-FPPZOLD ) 

40  FPPZOLD=FPPZNEW 
FP IOLD=FP I  NEW 
GO  TO  10 

****■»-»*■»*  ■*•*-*  3 

*  SOLVE  FOR  THE  DISTRIBUTION  OF  R  * 


100  L  =0 

C 1 =- 2  0*PP/< < UWUE-1  )**2) 


E2 


INITIALIZE  DATA  FOR  IMSL  SUBROUTINE  ' DGE AR ' 


N-5 

T  OL=  OOOOOl 
MZTH=2 
MI  TER  -3 

SET  THE  INITIAL.  CONDITIONS 

R Z  IS  R  '  0 > 

R  Z  =- 0 .  000825* CUWUE**2 )+0.  01155*UW UE+O  S799 
0  X  =0  0 

Y  (  1  '  =  0  0 

Y i2>  =  UWU£ 

V ( 3 ) =  FPP  Z 
Y < 4)-RZ 

Y  (  5)  =0.  0 
H  =  .  OOOOOl 
I NDEX  =  1 

SOLVE  FOR  K  AT  INFINITY  ( R I ) 

DO  120  1  =  1.-300 
X  END=  02*1 

CALL  DGEAR  (N,  FCN,  FCNJ.  X,  H,  Y,  XEND,  TOL,  METH,  MITER.  INDEX  ,  IWK,  WK,  IER 
STORE  THE  NON-DI MENS IONAL  BOUNDARY  LAYER  PARAMETERS  FOR  LATER  USE 
BLP(  1+1,  1  )=Y(  1  ) 

BLP  (  1  +  1,  2  >=Y<2) 

BLP  (  1+1,  3  >  =Y  (  3 > 

BLP  (1  +  1,  4  )=Y(4> 

BLP<  1+1,  5)=Y<  5) 

0  IF  ( IER  NE  0)  WR  I  TE  <  6,  *  )  '  IER2  IS  'TIER 

RZNEW=RZ 
R INEW=Y ( 4 ) 

I F ( IOUT  EQ  2)  WR I TE ( 6 , * )  'RINEW  IS  ' , R INEW, RZNEW 
L  =L+  1 

I F ( L  GT  15)  STOP 
I F ( L  GT  1 )  GO  TO  130 
R  Z=R  Z  NEW  +  O  05 
GO  TO  140 

CALCULATE  A  NEW  R  AT  ZERO  (RZ) 

iO  IF'ABSLRINEW-0  0)  LT  i.OE-05)  GO  TO  200 

R  Z  =R  Z  NEW+ (0  O-RINEW) / ( R I NEW-R I  OLD ) * ( R ZNEW-R ZOLD ) 

0  R  Z OLD  =P ZNEW 
R I  OLD  =R INEW 
GO  TO  110 


*  a-  a-a-a  -a-  a-  >-.i  -a-***  -#■  -a-  a- *  *  »•«■■»■  -a -a*-*  *■  * -turn-* 

*  SOLVE  FOR  THE  DISTRIBUTION  OF  S  * 
«■*•***  a-»a-a  -a  »  a -a  •#■•»■**■**•»•#■**■*■«■■#•»**•*■«■#■#••** 


200  R  EC=R  2 
L  =0 

C  1  =0  0 

INITIALIZE  DATA  FOR  IMSL  SUBROUTINE  'DGEAR ' 


N  =  5 

TOL=  000001 
METH  =  2 
M I  TER  =3 

SET  THE  INITIAL  CONDITIONS 
SPZ  IS  S  <  0  ■ 

SPZ  =  0  006  725*  OJWUE**2 >-0.  1 725 5*UWUE-0  533 
210  X=0  0 

Y  (  1  >=0.  0 
Y ( 2 ) =  UWUE 

Y  <  3) =FFFZ 

Y  <  4  )  =  1  0 
Y (5*=SPZ 
H=  000001 
I  NDEX  =  1 

SOLVE  FOR  S  AT  INFINITY  (SI) 

DO  220  1=1, 300 
X  END=  02*1 

CALL  DGEAR  ( N,  FCN,  FCNJ,  X,  H,  Y,  XEND,  TOL,  METH,  MITER,  INDEX,  INK,  WK,  IER  ) 
BLP  (  1*1,  fa  )  =Y  (  4  ) 

BLP(  1*1,  7)=Y(  5> 

220  IFUER  NE  0)  NR  I TE  ( 6,  *  )  '  IER3  IS  ',IER 

SP  ZNEU=SP Z 
S INEW=Y( 4 ) 

I F  ( I  OUT  EQ  2)  WR  I  TE  (  6 ,  * )  'SINEW  IS  '  ,  SINEW,  SIOLD,  SPZ 
L  =L+ 1 

I F  ( L  GT  15)  STOP 
I F (L.  GT  1  )  GO  TO  230 
SP Z  =  SP  ZNE W*0  05 
GO  TO  240 

CALCULATE  A  NEW  S-PR IME  AT  ZERO  (SPZ) 

230  IF  ( A  B  S  (  S I  NEW  -  0  0)  LT.  1.0E-05)  GO  TO  300 

SP  Z=SP  ZNE W* ( 0  O-SINEW ) / ( S INEW-S I OLD ) * ( SP ZNEW-SPZOLD ) 

240  SP  ZOL  D=3P  ZNEW 
S I  OLD  =SI NEW 
GO  TO  210 
300  CONTINUE 


E4 


LOAD  THE  BLF  MATRIX  FOR  THE  ROW  CORRESPONDING  TO  ETA 
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C  OUTPUT  THE  MATRIa  OF  NON -D I  MENS  I  ON  BOUNDARY  LAYER  PARAMETER 


DO  310  L-L301  5 
ETA-  02**  I -1) 

WRITE  -  a.  1000)  ETA,  (BLP(  I,  J),  J  =  l,  7  > 
310  CONTINUE 


C 

C  **********  !»■■■>>»  ********  ************  ****  ********  *  *** 
C  *  USE  BOUNDARY  LAYER  PROFILE  TO  DO  CALCULATIONS  * 
C  ****  ****  ************************  **** 


C  INITIALIZE  SOME  OF  THE  BOUNDARY  LAYER  PARAMETERS 

C 

320  RDN=0  C 
S  DN=0  0 

V  OL.D=  0  0 

C  1  = 1  G  AIR-i  )/  2.  *<  t  UWUE-1  )*ME)**2 
TWI=TE*(1  +C 1 *RtC  ) 

c 

C  REFERENCE  TO  A  MEAN  TEMPERATURE  TO  GET  A  BETTER  SOLUTION 

C  SEE  NACA  TN  3712 
C 

T  R  OF  =  0  5  »  i  TW+  TE )  •+  0  22*<TWI-TE  ) 

R  HOR  E F =R H  OE*T E / TR EF 
NUREF=VlSCO( TREF )  /RHOREF 
C 

C  CALCULATE  THE  BOUNDARY  LAYER  DATA 

C 

DO  40  0  -2,  30  1 

E  T  A=  02  * i K-l  ) 

FP-BL  P  •  *. .  2  ) 

R -BLP  (K.  4  ) 

5=BLT  iK.  t>  ) 

RDN=RDN+  (  BLP  <  K,  4  )  +  BLP  (K-l  ,  4)  )  /2.  *  02 
SDN -SDN*-  i  BLP  (  K ,  6  )  +BLP  (K-l  ,  6)  )  /2.  *.  02 

C 

T  =  T£*  '  1  * C  1*R  +  < TW-TW I  )/TE*S) 

RHO=RHOE*  TE/ T 

C  2-TE  ThEF* ( ETA+C 1 *RDN+ < T W-TW I )/TE*SDN) 

Y  NEW=  SOR  T  <  2  * NUR E  F / UE ) *C 2 
D  Y  =YNEw- Y  OLD 


U=UE*FP 

C3=( U-UW ) /(UE-UW) 
IF(C3  LT  99)  N=K 


WHILE  CALCULATING  Y  AT  CHAR.  LENGTH 


DO  410  K=l,  301, 5 

410  WRITE  •  6,  3000)  A*  K .  1  ) ,  A  <K,  2  )  *1 000  ,  A(K,  3) ,  A(K,  4) , 
1 A  ( K  2  >  /DELTA.  A  <  K ,  5),  AIK,  6) 


RETURN 


C  FORMAT  STATEMENTS 

1000  FORMAT*  I  H  ,  F4  1.7(5*.  F7  4  )  • 

2000  FORMAT  *  1H 1 ,  •'* -^BOUNDARY  LAYER  PROFILE  DATA*** - VERSION  3  ,  ) 

2100  FORMAT  *  .  H  ,  UE  =  '  ,  F7  1,'  M/SEC',  T42,  '  TE  =  F7  1.  '  DEG  K  •  ) 

2200  FORMAT*  1M  ,  'Uw  -  ',F7  1,  '  M/SEC;  T42,  •  TW  =  ',F7  1.  '  DEG  K'  > 

2300  FORMAT  <  1H  ,  'UW/UL  -  F4  2.  T40,  '  TREF  =  '  ,  F7  1,  '  DEG  K'  > 

2400  FORMA  T ■  1 H  , 'UE+AE  =  '  , F7  1,'  M/SEC  '  > 

2500  FORMAT*  .AO,  'REYNOLDS  NUMBER  ( L=1  MM  )  -  ',E9  4) 

2550  FORMAT* 1H0,  DISPLACEMENT  THICKNESS  ( L= 1  MM )  -  '  ,  F5  4,’  MM') 
2600  FORMAT* 1H  ,  DISPLACEMENT  WEDGE  ANGLE  ( L= 1  MM )  =  ' ,  F6.  3,  '  DEG') 
2700  FORMAT*,  1H  ,  '  ETA  Y  U 


real  FUNCTION  UISCO(TEMP) 


C  THIS  CALCULATION  USES  EQUATION  68  FROM  NASA  TECHNICAL  REPORT 
C  TR  R -50 

C 

VISCO^l  462E— G6#SQRT ( TEMP ) / ( 1 .  +112.  /TEMP) 

C 

RETURN 

END 


